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ABSTRACT 

The  setting  of  alarm  levels  plays  a  vital  role  in  a  machinery  condition  monitoring  and 
diagnostic  system.  In  this  research,  two  approaches  to  setting  vibration  alarm  levels  using 
vibration  signals  produced  by  fire  pumps  are  presented  in  the  time  and  frequency  domains. 
In  the  time  domain,  the  cross  peak  analysis  (CPA)  is  proposed  to  extract  the  dominate  peak 
points.  The  distribution  of  these  cross  peak  points  is  found  to  have  a  lognormal  distribution 
and  can  be  normalized  to  a  Normal  distribution  in  the  VdB  domain.  The  computed  (i+2o 
value  in  the  VdB  domain  is  suggested  for  use  as  the  alarm  level.  In  the  frequency  domain, 
the  1/1  octave  band  analysis  (OB  A)  is  introduced.  Three  artificial  fault  simulations  are  con- 
ducted to  compare  the  1/1  octave  band  method  with  the  broadband  method.  The  results 
show  that  the  1/1  octave  band  method  is  more  sensitive  to  the  changes  in  VdB  level  than 
the  broadband  method.  The  computer  programs  to  perform  these  two  analyses  are  written 
using  MATLAB.  Examples  of  the  use  of  these  programs  are  included  in  this  report. 
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I.  INTRODUCTION 

In  recent  years,  the  rising  cost  of  machine  maintenance  has  driven  engineers  to 
develop  more  economical  and  efficient  methods  to  determine  machine  health  and  to  plan 
accordingly  the  required  preventive  and  corrective  maintenance.  The  most  popular 
technique  in  use  today  is  a  predictive  maintenance  program  based  on  condition  monitoring. 
In  naval  applications,  condition  monitoring  is  commonly  achieved  utilizing  vibration 
measurement  and  analysis  on-board  surface  ships  and  submarines  [Refs.  1-7] 

This  research  focuses  on  the  use  of  vibration  measurement  to  monitor  machine  health 
and  to  diagnose  system  problems  which  could  lead  to  machine  failure.  In  addition  to 
providing  accurate  and  understandable  data  on  the  machine's  current  condition,  a 
monitoring  and  diagnostic  system  must  also  limit  the  number  of  false  alarms.  Alarm 
threshold  settings  are  vitally  important  in  machine  vibration  diagnostics.  Alarm  thresholds 
set  too  high  may  result  in  premature  machine  failure  caused  by  an  undetected  failure 
condition.  Alarm  thresholds  set  too  low  may  result  in  frequent  false  alarms  causing 
unnecessary  system  interruptions  and  repairs.  False  alarms  also  reduce  operator  confidence 
in  the  monitoring  and  diagnostic  system. 

Because  the  optimum  setting  of  alarm  thresholds  in  vibration  monitoring  and 
diagnostic  systems  continues  to  be  problematical,  the  goals  of  this  research  were: 

•  To  establish  a  statistical  analysis  method  for  setting  vibration  alarm  levels  by  using 
time  domain  data. 

•  To  investigate  the  difference  in  broadband  alarm  levels  between  the  cross  peak 
method  and  vibrometer  readings. 

•  To  perform  a  1/1  octave  band  analysis  for  setting  vibration  alarm  levels  by  using 
time  and  frequency  domain  data. 

Chapter  II  provides  the  background  information  on  the  machinery  utilized  in  this 
research  and  a  discussion  of  alarm  level  setting.  Chapter  III  describes  the  procedure  of  data 


acquisition  and  processing  and  the  methods  of  analysis.  In  the  time  domain,  a  statistical 
analysis  method  using  "cross  peak"  data  is  introduced.  In  the  frequency  domain,  a  1/1 
octave  band  concept  is  used  to  set  alarm  levels.  The  test  results  and  observations  using  time 
domain  data  are  presented  in  Chapter  IV.  Three  types  of  shipboard  fire  pumps  were 
analyzed  using  time  and  frequency  data.  Three  artificial  fault  simulations  were  performed 
to  assess  the  1/1  octave  band  method.  Chapter  V  demonstrates  the  1/1  octave  band  analysis 
results  using  the  NAVSSES  frequency  domain  database.  Chapter  VI  contains  conclusions 
and  recommendations  for  further  research. 


II.  BACKGROUND 

A.      DATA  SOURCE 

The  fire  pump  on  a  naval  ship  provides  more  than  fire  fighting  water  to  damage  control 
systems.  The  seawater  provided  by  fire  pumps  is  also  used  by  vital  air  conditioning  and 
chill  water  systems,  main  drainage  systems  and  decontamination  sprinkler  systems.  At  least 
one  fire  pump  is  running  whenever  a  ship  is  underway  or  at  anchor.  In  this  research,  23  fire 
pumps  from  four  classes  of  ships  are  grouped  into  three  types  in  Table  1 . 

TABLE  1:  CLASSIFICATION  OF  FIRE  PUMP  DATA 


Type 

Date  of 
Survey 

Number  of 
Fire  Pumps 

Number  of 
Pickup  Locations 

Number  of 
Data  Sets 

I 

10/30/89 

6 

5 

30 

2/26/90 

6 

5 

30 

n 

N/A 

3 

5 

15 

ni 

11/24/89 

4 

5 

20 

3/10/89 

4 

5 

20 

The  vibration  velocity  signals  were  measured  by  means  of  transducers  strategically 
placed  at  pickup  locations  on  the  pumps.  Schematic  layouts  of  each  fire  pump  and  pickup 
locations  are  depicted  in  Figures  1  through  3.  A  list  of  abbreviation  for  pickup  locations  is 
shown  in  Table  2.  The  transducer  pickup  placement  is  uniaxial  with  one  radial  pickup  at 
each  bearing  and  one  axial  pickup  at  the  thrust  bearing.  This  arrangement  is  listed  in  the 
Vibration  Test  and  Analysis  Guide  (VTAG)  [Ref.8].  This  guide  provides  the  most  current 
technical  information  for  each  machine  in  a  particular  ship  class  and  includes  information 
on  pickup  locations,  operating  conditions  and  a  table  of  vibration  source  components. 


Table  3  through  Table  5  summarize  the  vibration  source  components.  These  tables  identify 
exciting  elements  within  each  machine  and  list  the  vibration  frequencies  generated  by  each 
element.  The  vibration  frequencies  were  normalized  as  multiples  of  the  machine's  rotation 
rate,  called  "orders".  The  time  waveforms  are  recorded  on  magnetic  tape  by  a  frequency 
modulated  (FM)  recorder.  Each  pickup  location  on  the  pumps  was  recorded  for  a  one 
minute  time  series  record.  Some  of  the  important  specifications  of  the  tape  recorder  are 
provided  in  Table  6. 

In  addition  to  the  data  tapes,  six  survey  dates  of  data  in  ASCII  form  were  received  from 
the  NAVSSES  database.  These  data  contain  the  frequency  spectrum  levels  measured  by 
vibrometer. 


TABLE  2:  ABBREVIATIONS  OF  PICKUP  LOCATION 


Abbreviation 

Description 

MB(FE) 

Motor  Bearing  (Free  End) 

MB(CE) 

Motor  Bearing  (Coupling  End) 

MB(CE/A) 

Motor  Bearing  (Coupling  End/ Axial  direction) 

PB(FE) 

Pump  Bearing  (Free  End) 

PB(CE) 

Pump  Bearing  (Coupling  End) 

PB(CE/A) 

Pump  Bearing  (Coupling  End/ Axial  direction) 

UMB 

Upper  Motor  Bearing 

LMB 

Lower  Motor  Bearing 

LMB(A) 

Lower  Motor  Bearing  (Axial  direction) 

UPB 

Upper  Pump  Bearing 

LPB 

Lower  Pump  Bearing 

MB(FE)  MB(CE)  PB(CE)  PB(FE) 


PB(FE/A) 


Figure  1  Schematic  Layout  of  a  Type  I  Fire  Pump 


TABLE  3:  VIBRATION  SOURCE  COMPONENTS  OF  TYPE  I  FIRE  PUMP 


Driver  (Motor) 

Driven  (Pump) 

Description 

Element 

Order 

Description 

Element 

Order 

Motor  Shaft  (Ref.) 

1 

Pump  Shaft 

1 

Fan  Blading 

5 

5 

Impeller  Vanes 

6 

6 

Slots 

54 

54 

Bearing 

FAGWT 

Bars 

44 

44 

Poles 

2 

2 

Bearing 

MRC310 

Bearing 

MRC311 

Operation  Condition  of  Motor  is  3595  RPM  (=59.9  Hz) 

Figure  2  Schematic  Layout  of  a  Type  II  Fire  Pump 


TABLE  4:  VIBRATION  SOURCE  COMPONENTS  OF  TYPE  II  FIRE  PUMP 


Driver  (Motor) 

Driven  (Pump) 

Description 

Element 

Order 

Description 

Element 

Order 

Motor  Shaft  (Ref.) 

1 

Pump  Shaft 

1 

Fan  Blading 

7 

7 

Impeller  Vanes 

5 

5 

Slots 

48 

48 

Bearing 

SKF306 

Bars 

38 

38 

Poles 

2 

2 

Bearing 

MRC310 

Bearing 

MRC311 

Operation  Condition  of  Motor  is  3555  RPM  (=59.3  Hz) 

Figure  3  Layout  of  a  Type  III  Fire  Pump 


TABLE  5:  VIBRATION  SOURCE  COMPONENTS  OF  TYPE  III  FIRE  PUMP 


Driver  (Motor) 

Driven  (Pump) 

Description 

Element 

Order 

Description 

Element 

Order 

Motor  Shaft  (Ref.) 

1 

Pump  Shaft 

1 

Poles 

2 

2 

Impeller  Vanes 

5 

5 

Bearing 

Ball  Bearing 

Bearing 

SKF  6307 

Operation  Condition  of  Motor  is  3570  RPM  (=59.5  Hz) 

TABLE  6:  SPECIFICATIONS  OF  FM  RECORDER 


Manufacture 

Dallas  Instruments  Inc. 

Model 

4800  FM  Recorder 

Frequency  Response 

2  to  5000  Hz  within  1  dB 

Signal  to  Noise  Ratio 

40  dB  rms 

Harmonic  Distortion 

No  harmonic  above  -40  dB 

Tape  Type 

TDK  AD  C60  cassette 

B.      Alarm  Level 

The  fundamental  question  that  must  be  answered  before  vibration  condition 
monitoring  can  be  used  as  a  diagnostic  tool  for  machinery  is  what  alarm  levels  identify 
operation-limiting  faults.  Generally,  the  vibrations  of  a  system  can  be  characterized  by  a 
reduced  data  set  in  various  domains.  The  criteria  used  to  set  alarm  levels  can  be  considered 
either  in  the  time  domain  or  the  frequency  domain,  as  discussed  below. 

1.     Time  Domain  Criteria 

a.     Vibration  Severity  Criterion  Method 

The  simplest  time  domain  method  is  the  vibration  severity  criterion.  The  root 
mean  square  (RMS)  value  of  vibration  velocity  is  usually  measured  and  compared  with 
vibration  severity  charts  [Ref.  9].  Figure  4  shows  the  vibration  severity  chart  that  was  first 
introduced  by  T.  C  Rathbone  in  1939  and  then  later  refined  by  the  Instrument  Research  and 
Development  Corporation.  Various  national  standard  organizations  have  also  published 
standards  forjudging  vibration  severity.  For  example,  International  Standards  Organization 
(ISO)  standards  2372  and  3945  provide  severity  guidelines  for  machinery. 
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Figure  4  General  Machinery  Vibration  Severity  Chart  [Ref.  9] 


The  disadvantages  in  using  these  criteria  are  [Ref.  10]: 

•  The  criteria  can  be  used  only  for  specific  types  of  machinery  at  a  standard  operating 
condition. 

•  No  diagnostic  information  is  provided. 

•  The  criteria  are  less  sensitive  during  the  early  stages  of  damage. 

Thus,  the  vibration  severity  criterion  method  can  be  only  used  as  a  rough 
indicator  of  machinery  health. 


b.     Amplitude  Probability  Criterion  Method 

A  More  sophisticated  method  was  developed  using  statistical  analysis  to 
examine  the  distribution  of  vibration  amplitudes.  The  amplitudes  used  can  be  either  peak- 
to-peak,  peak  or  RMS  readings  of  displacement,  velocity  or  acceleration.  Both  Campel 
[Ref.  1 1]  and  Murphy  [Ref.  12]  use  this  statistical  method  to  establish  alarm  levels  based 
on  the  mean  of  the  reading  plus  3  standard  deviations.  The  key  drawback  of  this  approach 
is  that  it  assumes  that  a  Normal  distribution  of  the  linear  readings  exists.  Vibration  readings 
have  a  "skewed"  rather  than  "Normal"  distribution  [Ref.  13]. 

2.     Frequency  Domain  Criteria 

a.  Broadband  Criterion  Method 

The  broadband  criterion  method  utilizes  a  vibrometer  which  can  add  all  the 
energy  dissipated  over  a  wide  frequency  range  (typically  10  to  10,000  Hz).  The  overall 
energy  is  typically  calculated  by  applying  the  RMS  summation  method  to  the  spectrum.  If 
the  overall  energy  level  exceeds  a  predetermined  level,  then  alarms  are  triggered.  This  is  a 
simple  and  effective  method  in  most  cases.  But  in  some  cases  may  be  inadequate  because 
it  is  less  sensitive  to  small  changes  associated  with  bearing  defects  in  the  presence  of  large 
amplitude  frequency  components.  This  limitation  can  be  overcome  by  using  the  octave 
band  or  narrow  band  method. 

b.  Octave  Band  Criterion  Method 

The  octave  band  criterion  method  is  often  used  in  acoustics  to  determine  the 
energy  level  changes  due  to  noise  and  vibration.  This  method  utilizes  a  constant  percentage 
bandwidth  to  divide  the  frequency  range  of  interest  into  several  bands  which  provide  more 
detailed  information  than  the  broadband  presentation  A  commonly  used  bandwidth  is  the 
one-third  octave  band.  Early  researchers  used  this  method  to  check  the  change  in  each  band 
level  and  determine  if  the  amplitudes  exceeded  normal  values  [Ref.  14-15]  Chapter  III 
discusses  this  method  in  detail 
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c.     Narrowband  Criterion  Method 

Since  the  broadband  and  the  octave  band  criteria  methods  lack  detailed 
vibration  information,  a  narrowband  criterion  method  is  gaining  popularity  The  bandwidth 
may  be  up  to  10%  of  the  center  frequency  range  (typically  10  to  1,000  Hz).  The  improved 
resolution  is  generally  up  to  400  or  800  lines  over  the  frequency  range  of  interest.  With 
knowledge  of  the  forcing  frequencies  of  rotating  components,  the  narrowband  data  is  very 
useful  for  diagnosing  specific  faults.  This  diagnosis  can  be  time  consuming;  it  is  also  often 
difficult  to  explain  the  source  of  some  of  the  peaks  in  the  spectra. 

Frequency  criteria  are  used  in  U.S.  Navy  surface  ship  to  establish  alarm 
settings  for  machinery  [Ref.  3].  Vibration  data  for  each  individual  machine  and  location  are 
processed  to  yield  mean  and  standard  deviation  values  from  frequency  spectra  data. 
Generally,  broadband  levels  are  used  as  a  screening  tools,  and  narrowband  levels  are  used 
as  diagnostic  tools. 
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III.  METHODS  OF  ANALYSIS 

A.      DATA  ACQUISITION  AND  PROCESSING  SYSTEM 

A  block  diagram  of  the  data  acquisition  and  analysis  system  used  is  depicted  in  Figure 
5.  The  Model  4800  FM  recorder  plays  back  the  machinery  vibration  data  tapes  to  generate 
an  analog  signal.  By  using  a  connector,  this  analog  signal  is  distributed  to  the  data 
acquisition  software  and  the  oscilloscope  simultaneously.  The  oscilloscope  controls  the 

quality  of  the  data  by  monitoring  the  signal  time  waveform.  The  EASYEST  LX  software, 
developed  by  Keithley  Asyst,  was  used  as  a  data  acquisition  system.  This  PC-based  data 
acquisition  system  allowed  us  to  digitize  analog  signals,  convert  the  units  from  voltage  to 
velocity  and  store  the  velocity  data  sequences  into  the  hard  disk  of  a  personal  computer. 
Finally,  MATLAB  software  was  used  to  retrieve  these  data  sequences  and  perform  data 
processing.  MATLAB  is  a  very  powerful  software  tool  for  signal  processing  All  programs 
used  in  this  research  were  written  using  MATLAB  The  synopsis,  on-line  help 
documentation  and  program  code  used  can  be  found  in  Appendix  A.  Examples  for  the  use 
of  these  programs  to  perform  the  analysis  are  included  in  Appendix  B. 

All  of  the  data  sets  were  sampled  at  a  sampling  frequency  of  10  kHz.  The  sampling 
duration  for  each  data  set  was  36.684  seconds 
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Figure  5  Block  Diagram  of  Data  Acquisition  And  Processing  System 


13 


B.      TIME  DOMAIN  -  CROSS  PEAK  ANALYSIS 

It  is  obvious  that  the  peak  envelope  distribution  of  vibration  time  domain  signals  is 
relevant  to  the  vibration  severity  of  machinery  component.  Vibration  severity  criteria 
which  uses  measured  RMS  values  to  determine  the  severity  of  component  vibration 
without  considering  the  dispersion  of  the  signal  is  truly  a  rough  guess.  A  better  approach  is 
to  represent  the  severity  of  component  vibration  in  terms  of  a  percentage  acceptance  level. 
The  percentage  acceptance  level  provides  the  percentage  of  outcomes  which  will  not 
exceed  this  level  threshold. 

For  a  Normal  distribution,  the  acceptance  level  is  closely  related  to  the  mean  (|i)  and 
standard  deviation  (a).  For  example,  fi+  1.96a  corresponds  to  a  97.5%  acceptance  level 
which  means  only  a  2.5%  probability  the  signal  amplitude  will  exceed  the  |i+l  .96a  value. 
For  the  sake  of  simplicity,  the  most  popular  assumption  to  assume  the  probability 
distribution  of  the  signal  was  a  Gaussian  (or  Normal)  distribution.  However,  for  rotational 
machinery  vibration,  since  we  are  only  concern  with  the  magnitude  of  the  signal,  the 
probability  density  function  has  a  skewed  rather  than  a  Normal  distribution.  At  this  point, 
it  is  worthwhile  to  review  the  statistical  moments  and  central  moments  of  a  random  variable 

[Ref.16].  For  the  discrete  random  variable  x  with  probability  Pr(x),  the  mh  moment  E(x") 
is  defined  as: 

£(*")  =  5>*  Pr(xk)  (J) 

k=  1 

The  first  moment  is  particularly  important  in  many  applications  and  is  given  the  name  mean 
value  (ji): 

k=  1 

Of  greater  significance  are  the  central  moments  which  defined  as: 
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£[(*-u,)n]    =  ^(xk-nxy  Pr{xk)  (3) 


k=  1 


The  central  moment  for  n=  1  is  zero.  The  central  moment  for  n=2  is  a  very  important  which 
quantity  called  the  variance  of  the  random  variable; 

00 

E[(x-»x)2]    =  £(*,-ux)2  Pr{xk)  (4) 

k=  1 

The  standard  deviation,  which  corresponds  to  the  dispersion  of  the  random  variable,  can  be 
obtained  by  taking  the  square  root  of  the  variance,i.e., 


a  =  J  Variance  (5) 

As  noted  earlier,  the  mean  and  standard  deviation  are  often  used  to  compute  the  acceptance 
level. 

The  higher  order  central  moments  (n>2)  are  often  normalized  by  dividing  by  the  nth 
power  of  the  standard  deviation.  The  third  and  fourth  normalized  central  moments  are 
primarily  used  to  indicate  the  shape  of  the  probability  density  function.  They  are  defined 
by  the  following  equations: 

£[(*-ux)3] 

(6) 


Skewness  = 

3 

a 

Kurtosis  = 

£[(x-u 

4 

a 

y] 

(7) 

The  skewness  shows  information  about  the  position  of  the  peak  density  relative  to  the  mean 

value.  The  kurtosis  indicate  the  spread  in  the  distribution.  For  a  perfect  Normal  distribution, 

the  skewness  is  zero  and  the  kurtosis  is  three. 

Since  we  are  interested  in  the  peak  envelope  distribution  of  the  vibration  signal,  we 
use  the  cross  peak  data  points  instead  of  the  overall  peak  data  points.  The  cross  peak 
technique  can  extract  the  dominant  component  from  a  complicated  time  waveform  data  set 
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by  selecting  the  maximum  peak  value  between  every  zero  crossing.  Figure  6  illustrates  the 
difference  between  overall  peak  and  cross  peak  points. 

As  we  will  see  in  Chapter  IV,  the  probability  distributions  of  cross  peak  points  have 
highly  positive  skewness  and  can  be  approximated  as  a  lognormal  distribution.  This  means 
that  the  distributions  will  become  normal  if  the  data  is  transformed  by  a  log  or  natural  log 
algorithm.  The  method  chosen  to  accomplish  this  transformation  was  to  convert  the  linear 
velocity  reading  to  velocity  dB  (VdB)  reading  The  VdB  is  defined  as: 

Velocity  (VdB)   =  20  log—  (8) 

Vref 

Q 

where  Vre/\s  the  reference  level,  normally  0VdB=10    m/sec. 

Figure  7  shows  a  flow  chart  for  our  cross  peak  analysis  technique.  First,  the  sampled 
data  is  imported  into  a  MATLAB  worksheet.  The  DC  offset  is  then  removed  by  subtracting 
the  mean  value  of  the  sampled  data.  The  cross  peak  points  between  every  zero  crossing  is 
then  found  using  a  subroutine.  As  can  be  seen  from  Table  1,  the  dynamic  range  of  the  tape 
recorder  is  40  dB  This  means  that  amplitude  smaller  than  1%  of  the  maximum  peak  may 
be  contaminated  by  noise.  For  this  reason,  we  set  a  1%  threshold  to  eliminate  these 
contaminated  points.  Equation  8  was  then  used  to  transform  the  linear  velocity  scale  into 

the  VdB  domain.  The  reference  OVdB  was  10  m/sec.  After  we  computed  the  mean, 
standard  deviation,  skewness  and  kurtosis,  plots  of  the  probability  density  function  were 
generated  Finally,  the  statistical  analysis  report  can  be  sent  to  laser  printer. 
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Figure  7  Cross  Peak  Analysis  Flow  Chart 
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C.     FREQUENCY  DOMAIN  -  1/1  OCTAVE  BAND  ANALYSIS 

This  method  divides  the  frequency  spectrum  into  constant  percentage  bands  having 
the  same  ratio  of  bandwidth  to  center  frequency  [Ref.  17].  Each  band  has  an  upper 
frequency  limit  (/})  and  lower  frequency  limit  (/}).  The  center  frequency  (fc)  of  any  band  is 
defined  as 


f  =  Jffi  (V 

The  ratio  of  center  frequencies  of  successive  proportional  bands  is  the  same  as/?//}  for 

any  one  band  i.e., 

f     f 


=  r  -  (io) 

f     f     if 

One  third  of  an  octave  is  defined  as  n=  1/3.  The  American  National  Standards  Institute 
(ANSI)  preferred  center  frequencies  and  pass  bands  for  1/3  octave  bands  is  tabulated  at 
Table  7. 

The  method  chosen  to  examine  the  frequency  domain  alarm  levels  in  this  research  was 
to  use  1/1  octave  bands  (n=\).  For  the  10-5,000  Hz  frequency  range  of  our  system,  the 
frequency  spectrum  can  be  divided  into  9  bands  using  the  ANSI  preferred  center 
frequencies.  The  center  frequencies  and  pass  bands  covering  the  frequency  range  10-5,000 
Hz  in  1/1  octave  bands  are  given  in  Table  8. 
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TABLE  7:  ANSI  PREFERRED  CENTER  FREQUENCIES  AND  PASS  BANDS 

FOR  1/3  OCTAVE  BAND 


Center 
Frequency  (Hz) 

Pass  Band  (Hz) 

Lower  Limit 

Upper  Limit 

12.5 

12 

14 

16 

14 

18 

20 

18 

22.4 

25 

22.4 

28 

31.5 

28 

35.5 

40 

35.5 

45 

50 

45 

56 

63 

56 

71 

80 

71 

90 

100 

90 

112 

125 

112 

140 

160 

140 

180 

200 

180 

224 

250 

224 

280 

315 

280 

355 

400 

355 

450 

500 

450 

560 

630 

560 

710 

800 

710 

890 

1000 

890 

1120 

1250 

1120 

1400 

1600 

1400 

1800 

2000 

1800 

2240 

2500 

2240 

2800 

3150 

2800 

3550 

4000 

3550 

4500 

5000 

4500 

5600 

6300 

5600 

7100 

8000 

7100 

9000 

10000 

9000 

11200 
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TABLE  8:  ANSI  PREFERRED  CENTER  FREQUENCIES  AND  PASS  BANDS 

FOR  1/1  OCTAVE  BAND 


Center 
Frequency  (Hz) 

Pass  Band  (Hz) 

Lower  Limit 

Upper  Limit 

16 

11.2 

22.4 

31.5 

22.4 

45 

63 

45 

90 

125 

90 

180 

250 

180 

355 

500 

355 

710 

1000 

710 

1400 

2000 

1400 

2800 

4000 

2800 

5600 

8000 

5600 

11200 

Figure  8  shows  the  flow  chart  for  our  1/1  octave  band  analysis  using  the  MATLAB 
program  developed.  This  technique  used  frequency  domain  data  to  compute  the  1/1  octave 
band  levels.  In  order  to  have  better  resolution  at  lower  frequency,  it  was  necessary  to  use 
two  spectrum  data  with  different  sampling  frequencies.  For  the  time  domain  data,  the  lower 
sampling  frequency  spanned  2  kHz  and  the  higher  sampling  frequency  spanned  10  kHz. 
Twenty  blocks  of  data  is  sampled  at  each  sampling  frequency,  a  block  of  data  being  2048 
data  points.  After  sampling,  a  FFT  was  performed  for  each  block  with  a  Hanning  window 
to  obtain  the  smoothed  linear  velocity  spectrum.  The  averaging  VdB  spectrum  was  then 
produced  by  transforming  the  linear  spectrum  to  the  VdB  domain  and  taking  an  average 
over  the  20  data  blocks.  For  the  frequency  domain  data  received  from  the  NAVSSES 
database,  10  order  and  100  order  data  sets  were  used  as  the  lower  frequency  spectrum  and 
higher  frequency  spectrum,  respectively.  The  1st  through  6th  1/1  octave  band  levels  were 
then  computed  using  the  lower  frequency  sampling  spectrum  and  the  7th  through  9th  1/1 
octave  band  levels  obtained  using  the  higher  frequency  sampling  spectrum.  Finally,  the  1/ 
1  octave  band  analysis  was  done  after  combining  these  two  ranges  into  9  band  levels. 
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Figure  8  Flow  Chart  of  1/1  Octave  Band  Analysis  Technique 
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IV.  RESULTS  OF  USING  TIME  WAVEFORM  TAPES 

A.     CROSS  PEAK  ANALYSIS 

1.     Probability  Distribution  of  Cross  Peaks 

The  available  data  was  digitized  and  then  grouped  into  three  types  after 
examination  of  the  spectrum  patterns.  Typical  time  series  and  frequency  spectrums  for  each 
type  of  fire  pump  are  depicted  in  Figure  9  through  Figure  1 1 .  The  data  collection  sheets  are 
summarized  in  Table  9  through  Table  13  for  each  class  of  ship.  One  set  of  data  is  missing; 
the  actual  number  of  data  sets  analyzed  was  134.  The  operation  speed  of  the  shaft  was 
selected  as  a  reference  speed.  This  speed  was  used  to  normalize  the  frequency  scale  into 
"orders".  The  "B.B.  Level"  is  the  broadband  level  measured  by  vibrometer  when  the  data 
was  recorded.  This  level  represents  the  total  energy  of  the  velocity  reading  for  each 
measurement.  The  full  scale  represent  the  calibration  factor.  The  "B.B.  Aim  Level"  is  the 
broadband  alarm  level  set  by  NAVSSES  using  the  vibrometer  reading  data.  Both  the  "B.B. 
Level"  and  "B.B.  Aim  Level"  was  measured  for  a  frequency  range  from  10  Hz  to  10,000 
Hz.  The  remarks  column  describes  the  quality  of  the  data  set. 

In  order  to  understand  the  distributions  of  the  cross  peak  points,  probability 
histograms  for  each  set  of  data  were  plotted  using  linear  velocity  scale.  Figure  12  and 
Figure  13  give  an  example  of  probability  histograms  for  type  I  fire  pumps  measured  at 
MB(FE).  Inspection  of  these  histograms  shows  the  distributions  to  have  an  exponential 
shape  for  the  most  part.  With  a  linear  velocity  scale,  this  type  of  distribution  means  that  the 
population  is  mainly  dominated  by  small  peaks  which  is  an  expected  result  for  machinery 
in  good  condition.  Those  distributions  having  a  large  population  at  higher  amplitudes  are 
regarded  as  representing  damaged  pumps.  Examination  of  these  damaged  data  sets  shows 
their  broadband  levels  to  be  higher  than  the  others  and  closer  to  the  broadband  alarm  levels. 
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Figure  9  Typical  Time  Series  and  Frequency  Spectrum  of  Type  I  Fire  Pumps 

(10/30/89)  Measured  at  MB(FE) 
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Figure  10  Typical  Time  Series  and  Frequency  Spectrum  of  Type  II  Fire  Pumps 

Measured  at  MB(FE) 
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Figure  1 1  Typical  Time  Series  and  Frequency  Spectrum  of  Type  III  Fire  Pumps 

(3/10/89)  Measured  at  UMB 
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TABLE  9:  DATA  COLLECTION  SHEET  FOR  TYPE  I  (10/30/89)  FIRE  PUMPS 

Machine 
Name 

Speed 
(RPM) 

Pickup 
Location 

B.B. 

Level 

(VdB) 

Full 
Scale 
(VdB) 

B.B.  Aim 
Level 
(VdB) 

Remarks 

Fire  Pump 

#1 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

102 

120 

109 

PB(CE) 

110 

120 

118 

PB(FE) 

111 

120 

119 

PB(FE/A) 

114 

120 

121 

Fire  Pump 

#2 

3595 

MB(FE) 

97 

100 

107 

MB(CE) 

101 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

106 

120 

119 

PB(FE/A) 

107 

120 

121 

Fire  Pump 
#3 

3595 

MB(FE) 

97 

100 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

111 

120 

118 

PB(FE) 

109 

120 

119 

PB(FE/A) 

113 

120 

121 

Fire  Pump 

#4 

3595 

MB(FE) 

100 

120 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

110 

120 

119 

PB(FE/A) 

111 

120 

121 

Fire  Pump 
#5 

3595 

MB(FE) 

100 

120 

107 

MB(CE) 

101 

120 

109 

PB(CE) 

110 

120 

118 

PB(FE) 

111 

120 

119 

PB(FE/A) 

114 

120 

121 

Fire  Pump 
#6 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

113 

120 

118 

Damaged 

PB(FE) 

111 

120 

119 

PB(FE/A) 

111 

120 

121 
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TABLE  10:  1 

DATA  COLLECTION  SHEET  FOR  TYPE  I  (2/26/90)  FIRE  PUMPS 

Machine 

Name 

Speed 
(RPM) 

Pickup 
Location 

B.B. 

Level 

(VdB) 

Full 
Scale 
(VdB) 

BB  Aim 
Level 

(VdB) 

Remarks 

Fire  Pump 
#1 

3595 

MB(FE) 

98 

100 

107 

MB(CE) 

97 

100 

109 

PB(CE) 

106 

120 

118 

PB(FE) 

110 

120 

119 

PB(FFVA) 

111 

120 

121 

Fire  Pump 
#2 

3595 

MB(FE) 

97 

120 

107 

MB(CE) 

96 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

107 

120 

119 

PB(FE/A) 

108 

120 

121 

Fire  Pump 
#3 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

116 

120 

118 

Damaged 

PB(FE) 

115 

120 

119 

Damaged 

PB(FFVA) 

117 

120 

121 

Fire  Pump 

#4 

3595 

MB(FE) 

98 

120 

107 

MB(CE) 

97 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

111 

120 

119 

PB(FFVA) 

112 

120 

121 

Fire  Pump 
#5 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

101 

120 

109 

PB(CE) 

109 

120 

118 

PB(FE) 

110 

120 

119 

PB(FE/A) 

113 

120 

121 

Fire  Pump 
#6 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

98 

120 

109 

PB(CE) 

115 

120 

118 

Damaged 

PB(FE) 

112 

120 

119 

PB(FE/A) 

113 

120 

121 
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TABLE  11:  DATA  COLLECTION  SHEET  FOR  TYPE  II  FIRE  PUMPS 


Machine 
Name 

Speed 
(RPM) 

Pickup 
Location 

B.B. 

Level 

(VdB) 

Full 
Scale 
(VdB) 

BB  Aim 
Level 
(VdB) 

Remarks 

Fire  Pump 
#1 

3555 

MB(FE) 

97 

120 

108 

MB(CE) 

99 

120 

112 

Bad  Data 

MB(CE/A) 

110 

120 

112 

PB(CE) 

110 

120 

117 

PB(FE) 

109 

120 

119 

Fire  Pump 
#2 

3555 

MB(FE) 

91 

100 

108 

MB(CE) 

99 

120 

112 

MB(CE/A) 

107 

120 

112 

PB(CE) 

108 

120 

117 

PB(FE) 

109 

120 

119 

Fire  Pump 
#3 

3555 

MB(FE) 

102 

120 

108 

MB(CE) 

103 

120 

112 

MB(CE/A) 

112 

120 

112 

PB(CE) 

113 

120 

117 

PB(FE) 

116 

120 

119 
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TABLE  12:  DATA  COLLECTION  SHEET  FOR  TYPE  III  (11/24/89)  FIRE  PUMPS 

Machine 
Name 

Speed 
(RPM) 

Pickup 
Location 

B.B. 

Level 

(VdB) 

Full 
Scale 
(VdB) 

BB  Aim 
Level 
(VdB) 

Remarks 

Fire  Pump 
#1 

3570 

UMB 

110 

120 

108 

Damaged 

LMB 

111 

120 

110 

Damaged 

LMB(A) 

103 

120 

108 

UPB 

117 

120 

120 

Damaged 

LPB 

112 

120 

121 

Fire  Pump 
#2 

3570 

UMB 

97 

120 

108 

LMB 

98 

120 

110 

LMB(A) 

98 

120 

108 

UPB 

107 

120 

120 

LPB 

105 

120 

121 

Fire  Pump 
#3 

3570 

UMB 

101 

120 

108 

LMB 

104 

120 

110 

LMB(A) 

101 

120 

108 

UPB 

109 

120 

120 

LPB 

111 

120 

121 

Fire  Pump 
#4 

3570 

UMB 

98 

120 

108 

LMB 

103 

120 

110 

LMB(A) 

98 

120 

108 

UPB 

114 

120 

120 

LPB 

113 

120 

121 
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TABLE  13:  DATA  COLLECTION  SHEET  FOR  TYPE  III  (3/10/89)  FIRE  PUMPS 


Machine 
Name 

Speed 
(RPM) 

Pickup 
Location 

B.B. 
Level 
(VdB) 

Full 
Scale 
(VdB) 

BB  Aim 
Level 
(VdB) 

Remarks 

Fire  Pump 
#1 

3570 

UMB 

98 

120 

108 

LMB 

99 

120 

110 

LMB(A) 

97 

120 

108 

UPB 

110 

120 

120 

LPB 

106 

120 

121 

Fire  Pump 
#2 

3570 

UMB 

96 

120 

108 

LMB 

104 

120 

110 

Missing 

LMB(A) 

102 

120 

108 

UPB 

117 

140 

120 

LPB 

114 

140 

121 

Fire  Pump 
#3 

3570 

UMB 

105 

120 

108 

LMB 

102 

120 

110 

LMB(A) 

101 

120 

108 

UPB 

111 

120 

120 

LPB 

107 

120 

121 

Fire  Pump 
#4 

3570 

UMB 

96 

100 

108 

LMB 

102 

120 

110 

LMB(A) 

99 

120 

108 

UPB 

114 

120 

120 

LPB 

109 

120 

121 
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Probability  Density  Function  for  Cross  Peak  Points 

Type  I  Fire  Pump  (10/30/89) 

Pickup  Location  :  MB(FE) 
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Figure  12  The  Probability  Distributions  of  Cross  Peak  Data  for  Type  I  (10/30/89)  Fire 

Pumps  Measured  at  MB(FE) 
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Probability  Density  Function  for  Cross  Peak  Points 
Type  I  Fire  Pump  (2/26/90) 
Pickup  Location  :  MB(FE) 
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Figure  13  The  Probability  Distributions  of  Cross  Peak  Data  for  Type  I  (2/26/90)  Fire 

Pumps  Measured  at  MB(FE) 
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2.     Statistical  Analysis  Results 

As  mentioned  before,  the  distributions  of  the  cross  peak  points  in  a  linear  velocity 
scale  have  an  exponential  shape  histogram.  This  shape  histogram  can  be  approximated  as 
a  lognormal  distribution.  Because  a  velocity  decibel  is  a  unit  of  measurement  equal  to  the 
logarithm  of  the  ratio  V/Vrej,  it  is  convenient  to  use  the  VdB  transformation  to  convert  an 
exponential  shaped  distribution  into  a  Normal  distribution  Before  we  applied  this 
transform,  the  1%  threshold  was  used  to  remove  distorted  data.  The  data  sets  with  same 
type  of  fire  pump  and  same  pickup  location  were  then  added  together  to  form  a  combined 
data  set  to  perform  the  statistical  analysis.  The  results  are  collected  in  Appendix  C.  Figure 
C.l  through  Figure  C  15  present  the  statistical  analysis  results  for  all  of  the  available  data. 
Figure  C.16  through  Figure  C.21  show  the  statistical  analysis  results  without  including  the 
damaged  data  sets.  Table  14  tabulates  these  results  in  the  linear  velocity  and  VdB  domains. 
A  comparison  of  alarm  thresholds  between  the  broadband  alarm  levels  (B.B.  Aim.  Level) 
and  the  computed  \jl+2o  level  is  shown  in  Table  15.  It  should  be  mentioned  that  these  two 
alarm  thresholds  are  based  on  different  methods  using  data  from  different  domains.  The 
broadband  alarm  level,  as  noted  earlier,  is  obtained  by  using  vibrometer  readings  in  the 
frequency  domain.  The  computed  u.+a  levels  represent  the  97.5%  acceptance  level  of  the 
cross  peak  envelopes  in  the  VdB  domain.  Several  observations  can  be  made: 

•  The  distribution  of  the  cross  peak  points  in  the  linear  velocity  domain  have  a  highly 
positive  skew.  This  can  be  treated  as  a  lognormal  distribution. 

•  By  using  the  Vdb  transformation,  the  distributions  in  the  linear  velocity  domain  can 
be  transformed  to  Normal  or  near  Normal  distributions.  Generally  speaking,  the  data  sets 
at  pickup  locations  for  the  motor  have  a  skewness  less  than  0. 1  and  the  data  sets  at  pickup 
locations  for  the  pump  ends  have  a  skewness  varying  between  0. 1  and  0.5. 

•  As  can  be  seen  from  Table  14,  removing  the  damaged  data  sets  does  not  cause  the 
probability  distributions  in  the  Vdb  domain  to  become  much  closer  to  Normal  as  measured 
by  the  skewness  and  kurtosis  For  the  MB(CE)  pickup  location  of  a  type  II  pump,  not 
considering  the  damaged  pump  data  sets  makes  the  skewness  even  more  severe  than  when 
the  damaged  pump  data  sets  are  considered.  Since  there  were  only  3  type  II  pumps,  this 
result  is  not  surprising 

•  In  general,  the  B.B.  Alarm  Levels  and  computed  u+2o  levels  in  the  VdB  domain  are 
quite  close.  The  relative  errors  between  these  two  methods  for  Type  I  fire  pumps  are  less 
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than  5%.  It  is  also  noted  from  Table  15  that  the  computed  u.+2o  levels  in  cases  with  small 
skewness  is  closer  to  the  broadband  alarm  levels  than  in  cases  with  high  skewness.  This  is 
because  the  computed  \±+2o  levels  for  the  Normal  distribution  is  very  sensitive  to 
skewness. 
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TABLE  14:  STATISTICAL  ANALYSIS  RESULTS  FOR  CROSS  PEAK  POINTS 


Type 

Pickup 
Position 

Linear  Velocity  Domain 

VdB  Domain 

(mm/s) 

a 

(mm/s) 

Skew. 

Kurt. 

(VdB) 

a 

(VdB) 

Skew. 

Kurt. 

I 

MB(FE) 

0.495 

0.533 

1.850 

6.201 

89.36 

9.20 

-0.046 

2.378 

MB(CE) 

0.539 

0.608 

2.052 

7.195 

90.20 

9.12 

0.049 

2.444 

PB(CE) 

2.828 

2.680 

1.517 

5.638 

104.6 

9.78 

-0.444 

2.476 

PB(FE) 

2.717 

2.681 

1.253 

3.959 

103.6 

10.45 

-0.335 

2.232 

PB(FE/A) 

3.699 

3.307 

1.086 

3.515 

106.9 

9.90 

-0.519 

2.472 

PB(CE)* 

2.319 

1.988 

1.046 

3.304 

103.3 

9.25 

-0.510 

2.506 

PB(FE)* 

2.525 

2.472 

1.169 

3.409 

103.0 

10.30 

-0.317 

2.228 

II 

MB(FE) 

0.504 

0.604 

2.059 

6.976 

88.4 

10.50 

-0.139 

2.420 

MB(CE) 

0.727 

0.763 

3.243 

66.02 

93.1 

8.78 

-0.022 

2.320 

MB(CE/A) 

2.148 

2.493 

1.321 

3.604 

100.0 

11.49 

0.019 

1.926 

PB(CE) 

2.484 

2.344 

0.905 

2.675 

102.6 

10.84 

-0.368 

2.038 

PB(FE) 

2.166 

2.739 

2.008 

6.446 

100.8 

10.29 

0.113 

2.366 

MB(CE)* 

0.619 

0.811 

2.146 

6.762 

90.5 

9.19 

0.495 

2.566 

III 

UMB 

0.672 

0.926 

4.985 

36.52 

92.1 

8.69 

0.029 

2.953 

LMB 

0.923 

1.046 

3.350 

19.59 

95.3 

8.56 

-0.077 

2.697 

LMB(A) 

0.703 

0.648 

1.881 

7.538 

93.3 

846 

-0.349 

2.663 

UPB 

4.141 

4.442 

1.570 

5.446 

106.1 

11.84 

-0.408 

2.207 

LPB 

3.160 

3.246 

1.779 

6.775 

105.1 

10.20 

-0.373 

2.430 

UMB* 

0.597 

0.602 

2.488 

11.46 

91.8 

8.34 

-0.167 

2.595 

LMB* 

0.776 

0.736 

1.980 

7.430 

94.3 

8.15 

-0.181 

2.555 

UPB* 

3.603 

3.906 

1.711 

6.361 

105.1 

11.60 

-0.367 

2.176 

:  The  statistical  analysis  results  without  considering  the  damaged  data  sets. 
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TABLE  15:  COMPARISON  OF  COMPUTED  n+2a  LEVEL  AND  BROADBAND 

ALARM  LEVEL 


Type 

Pickup 
Position 

Skewness 

(VdB) 

a 

(VdB) 

H+2a 
(VdB) 

B.B.  Alarm 
Level 

(VdB) 

Relative 
Error 

I 

MB(FE) 

-0.046 

89.36 

9.20 

107.8 

107 

0.71% 

MB(CE) 

0.049 

90.20 

9.12 

108.3 

109 

-0.68% 

PB(CE) 

-0.444 

104.6 

9.78 

124.2 

118 

4.95% 

PB(FE) 

-0.335 

103.6 

10.45 

124.5 

119 

4.42% 

PB(FE/A) 

-0.519 

106.9 

9.90 

126.7 

121 

4.50% 

PB(CE)* 

-0.510 

103.3 

9.25 

121.8 

118 

3.11% 

PB(FE)* 

-0.317 

103.0 

10.30 

123.6 

119 

3.72% 

II 

MB(FE) 

-0.139 

88.4 

10.50 

109.4 

108 

1.27% 

MB(CE) 

-0.022 

93.1 

8.78 

110.6 

112 

-1.25% 

MB(CE/A) 

0.019 

100.0 

11.49 

123.0 

112 

8.91% 

PB(CE) 

-0.368 

102.6 

10.84 

124.3 

117 

5.86% 

PB(FE) 

0.113 

100.8 

10.29 

121.4 

119 

1.96% 

MB(CE)* 

0.495 

90.5 

9.19 

108.9 

118 

-8.37% 

in 

UMB 

0.029 

92.1 

8.69 

109.5 

108 

1.40% 

LMB 

-0.077 

95.3 

8.56 

112.4 

110 

2.11% 

LMB(A) 

-0.349 

93.3 

8.46 

110.2 

108 

2.04% 

UPB 

-0.408 

106.1 

11.84 

129.8 

120 

7.54% 

LPB 

-0.373 

105.1 

10.20 

125.5 

121 

3.59% 

UMB* 

-0.167 

91.8 

8.34 

108.5 

108 

0.44% 

LMB* 

-0.181 

94.3 

8.15 

110.6 

110 

0.52% 

UPB* 

-0.367 

105.1 

11.60 

128.3 

120 

6.47% 

:  The  statistical  analysis  results  without  considering  the  damaged  data  sets. 
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B.      1/1  OCTAVE  BAND  ANALYSIS 

1.     1/1  Octave  Band  Analysis  Results 

The  1/1  octave  band  analysis  was  performed  for  all  types  of  fire  pumps  with  the 
time  domain  data.  Band  1  through  band  6  were  computed  by  using  fs=2  kHz  to  provide 
better  resolution  at  lower  frequencies.  Band  7  to  band  9  were  computed  by  using  fs=10  kHz. 

Figure  14  through  Figure  16  show  the  summarized  1/1  octave  band  levels  for  each  type  of 
fire  pump  without  considering  the  damaged  data  sets.  The  10-5,000  Hz  broadband  levels 
are  also  shown  on  these  plots.  It  was  computed  by  adding  all  the  VdB  spectrum  using  a 
RMS  algorithm.  The  means  (|i)  and  standard  deviations  (a)  of  each  octave  band  levels  for 
the  same  type  fire  pumps  with  the  same  pickup  locations  was  also  computed.  The  thin  line 
represents  the  mean  (fi)  band  level  and  the  thick  line  represents  the  mean  plus  one  standard 
deviation  (uxj)  band  level.  Table  16  tabulates  these  band  levels. 

Generally  speaking,  the  dominate  levels  are  located  in  the  first  six  bands  (10  to 
710  Hz).  This  implies  that  the  energy  of  vibration  is  concentrated  at  lower  frequencies.  For 
fire  pumps  with  a  3555  RPM  operation  speed,  these  bands  are  approximated  up  to  12 
orders.  Upon  examination  of  the  third  band  through  the  sixth  band  for  each  data  set,  it  is 
obvious  that  the  band  levels  for  damaged  data  sets  are  higher  than  the  for  undamaged  pump 
data  sets. 
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1/1  Octave  Band  Analysis 

Type  I  Fire  Pump 
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Figure  14  The  Summarized  1/1  Octave  Band  Levels  for  Type  I  Fire  Pumps 
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Figure  14  The  Summarized  1/1  Octave  Band  Levels  for  Type  I  Fire  Pumps  (Cont.) 
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1/1  Octave  Band  Analysis 

Type  II  Fire  Pump 

Each  Data  Set  Use  20  Averaged  FFT  Data 
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Figure  15  The  Summarized  1/1  Octave  Band  Levels  for  Type  II  Fire  Pumps 
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(e)  Pickup  Location  :  PB(FE) 


Use  3  data  sets 
10-5k  Hz  B.B.  u-kt  LeveM  14.71  dB 

u:    

U.+C  :     
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Figure  15  The  Summarized  1/1  Octave  Band  Levels  for  Type  II  Fire  Pumps  (Cont) 
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1/1  Octave  Band  Analysis 

Type  III  Fire  Pump 

Each  Data  Set  Use  20  Averaged  FFT  Data 
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(b)  Pickup  Location  :  LMB 


Use  6  data  sets 
10-5k  Hz  B.B.  u-kj  Level=1 02.45  dB 
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Figure  16  The  Summarized  1/1  Octave  Band  Levels  for  Type  III  Fire  Pumps 
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(c)  Pickup  Location  :  LMB(A) 


Use  8  data  sets 
10-5k  Hz  B.B.  u+o  Level=101.73  dB 
u:     
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Frequency  (Hz) 

(d)  Pickup  Location  :  UPB 
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Use  7  data  sets 
10-5k  Hz  B.B.  h-kj  LeveM  15.25  dB 
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Figure  16  The  Summarized  1/1  Octave  Band  Levels  for  Type  III  Fire  Pumps  (Cont) 
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TABLE  16:  SUMMARIZED  MEAN  (n)  AND 

OCTAVE  BAND 


STANDARD  DEVIATION  (a)  OF  1/1 
LEVELS 


Type 

Pickup 

Band  Level  (VdB) 
0VdB=le-8  m/sec 

1st 

2nd 

3rd 

4th 

5th 

6th 

7th 

8th 

9th 

I 

MB(FE) 

H 

84.5 

75.5 

95.0 

88.3 

83.4 

81.3 

69.4 

70.3 

73.4 

a 

4.87 

3.59 

1.56 

1.09 

2.72 

3.34 

4.84 

4.16 

4.70 

MB(CE) 

H 

78.9 

81.8 

95.5 

91.3 

83.6 

83.0 

72.5 

71.2 

75.1 

o 

4.48 

5.49 

1.79 

2.54 

1.34 

2.42 

1.09 

1.95 

3.45 

PB(CE) 

H 

83.2 

85.9 

101.3 

98.4 

100.5 

99.0 

86.8 

83.4 

86.2 

a 

3.42 

5.05 

2.90 

4.12 

2.68 

3.29 

2.54 

4.11 

4.56 

PB(FE) 

[A 

85.0 

89.7 

102.9 

95.3 

103.9 

98.3 

84.0 

82.7 

85.3 

O 

4.76 

6.14 

2.90 

3.17 

3.60 

3.96 

3.71 

4.89 

3.38 

PB(FE/A) 

^ 

82.8 

86.3 

101.6 

99.1 

106.6 

102.3 

89.2 

88.9 

86.3 

o 

3.41 

3.21 

4.41 

5.30 

3.65 

4.29 

2.67 

3.86 

3.48 

II 

MB(FE) 

H 

74.0 

75.0 

90.0 

89.5 

83.7 

83.2 

69.8 

68.8 

71.7 

a 

7.67 

1.99 

12.44 

1.37 

2.39 

3.73 

5.28 

4.93 

8.23 

MB(CE) 

n 

73.6 

77.3 

98.1 

92.8 

85.2 

86.2 

72.2 

70.8 

76.5 

a 

6.03 

0.48 

2.14 

4.01 

0.12 

0.60 

0.05 

0.57 

0.94 

MB(CE/A) 

H 

81.3 

81.9 

102.9 

102.0 

89.8 

94.2 

82.0 

76.6 

80.2 

o 

1.86 

3.44 

11.06 

3.87 

1.71 

2.09 

3.09 

2.38 

1.72 

PB(CE) 

H 

86.4 

84.3 

106.3 

97.6 

91.2 

98.4 

83.9 

76.3 

80.3 

o 

1.36 

3.72 

6.00 

5.37 

1.18 

2.19 

2.19 

2.59 

1.93 

PB(FE) 

n 

78.2 

81.2 

108.6 

99.3 

92.3 

92.0 

85.9 

83.7 

82.1 

o 

3.55 

4.13 

5.57 

4.88 

0.95 

1.26 

1.38 

2.35 

3.18 

III 

UMB 

n 

87.6 

79.1 

91.1 

89.2 

89.4 

84.3 

73.9 

77.4 

76.1 

o 

3.67 

3.99 

5.29 

5.21 

3.70 

3.22 

4.96 

2.75 

6.76 

LMB 

n 

87.0 

81.7 

95.6 

94.0 

89.3 

85.7 

74.3 

79.7 

79.7 

o 

4.19 

3.38 

3.42 

2.95 

3.72 

2.49 

2.79 

1.46 

3.18 

LMB(A) 

n 

87.1 

79.3 

94.7 

91.9 

87.2 

85.9 

74.4 

75.3 

80.4 

a 

3.19 

3.35 

4.57 

2.71 

3.20 

2.93 

1.21 

4.06 

3.57 

UPB 

H 

91.4 

88.7 

106.3 

101.3 

99.3 

97.8 

85.2 

80.5 

83.9 

o 

4.55 

6.51 

6.57 

4.24 

8  18 

7.73 

7.86 

7.04 

7.98 

LPB 

H 

93.8 

93.3 

100.4 

100.6 

102.8 

99.1 

88.3 

83.2 

85.3 

o 

4.60 

4.66 

4.94 

2.21 

4.35 

4.83 

5.64 

6.99 

6.61 
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2.     Artificial  Fault  Simulations 

a.     Fault  Simulation 

In  order  to  assess  the  1/1  octave  band  method,  three  cases  of  artificial  fault 
simulation  were  performed.  Before  we  simulated  these  cases,  we  had  to  define  a  "fault". 
For  simplicity,  we  assumed  a  fault  could  be  approximated  with  5  bars,  as  shown  in  Figure 
17.  The  bandwidth  of  the  fault  is  about  4.5  Hz. 


Figure  17  Fault  with  6VdB  Gain  at  Center  Bar 

Based  on  Eq.  (8),  if  the  vibration  amplitude  doubles  with  a  linear  velocity 
scale  then  the  VdB  level  will  increase  6VdB  at  the  corresponding  frequency.  Therefore,  we 
defined  a  fault  as  having  a  6VdB  gain  at  the  center  bar  and  use  four  bars  with  3VdB  (half 
the  power  of  the  center  bar)  and  2  VdB  at  the  side  bands  to  simulate  leakage  effects  around 
the  corresponding  frequency. 
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b.     Simulation  of  A  Misalignment  Fault 

Misalignment  occurs  when  the  center  lines  of  two  shafts  are  offset  or  meet  at 
an  angle.  The  characteristics  of  misalignment  in  a  spectrum  include: 

•  High  amplitude  axial  peaks  and  radial  peaks  at  1,  2  and  3  orders  of  shaft  RPM. 

•  Higher  harmonics  of  the  shaft  RPM  (greater  than  4  orders  of  shaft  RPM)  are 
generally  low  in  amplitude. 

Thus,  it  is  reasonable  to  add  the  6VdB  gain  fault  at  the  first  and  second  order. 
Figure  1 8(a)  and  1 8(b)  show  the  spectrums  before  and  after  imposing  the  fault.  Figure  1 8(c) 
compares  the  differences  in  the  1/1  octave  band  levels  and  the  10-5,000  Hz  broad  band 
level  for  these  two  conditions.  As  can  be  seen  in  Figure  18(c),  the  1/1  octave  band  levels 
have  a  4.46  VdB  gain  at  third  band  (which  corresponds  to  1  order  of  shaft  RPM)  and  a  3.13 
VdB  gain  at  the  forth  band  (which  corresponds  to  2  orders  of  shaft  RPM)  However,  the  10- 
5,000  Hz  broadband  level  only  increased  2. 1 7  VdB.  Obviously,  the  1/1  octave  band  method 
is  more  sensitive  than  broadband  method. 
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Case  ^Misalignment  Fault 

Type  I  Fire  Pump 
Pickup  Location  :  PB(CE) 
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(c)  1/1  Octave  Band  Levels 
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Figure  18  Simulation  Result  for  Artificial  Misalignment  At  Shaft 
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c.     Simulation  of  A  Looseness  Fault  At  Impeller 

Impeller  looseness  is  a  rotating  element  looseness.  The  important 
characteristics  of  looseness  in  a  spectrum  include: 

•  Presence  of  a  large  number  of  harmonics. 

•  Presence  of  half-harmonics. 

For  example,  the  fire  pump  of  a  Type  I  has  six  impellers  contributing  to  a 
forcing  frequency  of  6  orders  of  the  shaft  RPM,  and  we  expect  higher  levels  at  6  orders  and 
its  harmonics  (12  orders,  24  orders,...  etc.)  for  a  loose  impeller.  The  levels  at  half- 
harmonics  (e.g.,  3  orders  and  9  orders)  will  also  increase.  Thus,  in  this  looseness  fault 
simulation,  an  artificial  fault  was  applied  by  adding  a  fault  with  6VdB  at  center  bar  at  the 
6th  order(i.e.,  6  x  shaft  RPM)  and  adding  two  fault  with  3  VdB  at  center  bar  at  the  3rd  order 
and  the  9th  order.  Figure  19(a)  and  19(b)  show  the  spectrum  of  good  and  damaged  fire 
pumps.  Figure  19(c)  compares  the  differences  in  1/1  octave  band  levels  and  the  10-5,000 
Hz  broadband  level  for  these  two  conditions.  At  the  fifth  of  the  1/1  octave  band  levels,  there 
is  a  1.5  VdB  gain.  However,  the  10-5,000  Hz  broadband  level  only  increased  0.29  VdB. 
Again,  the  1/1  octave  band  method  is  more  sensitive  than  broadband  method. 
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Case  II: Looseness  Fault  At  Impeller 

Type  I  Fire  Pump 
Pickup  Location  :  PB(FE) 


(a)  Good  Condition 


(b)  Damaged  Condition 
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(c)  1/1  Octave  Band  Levels 


10-5kHz  B.B.  Level 
—  Good   1 08.01  db 


102.08 

100.58 

90.84 

1 

90.81 

10  10 

Frequency  (Hz) 


10 


Figure  19  Simulation  Result  for  Artificial  Looseness  At  Impeller 
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d.    Simulation  of  A  Bearing  Fault 

For  a  steady  state  condition,  some  periodic  signatures  exist  which  relate  to 
corresponding  bearing  faults,  called  bearing  frequencies  These  bearing  frequencies  can  be 
found  in  the  VTAG.  In  this  section,  an  artificial  fault  imposed  at  the  motor  bearing  (free 
end)  of  a  Type  I  fire  pump  has  been  simulated.  The  bearing  frequencies  of  this  bearing  are 
tabulated  in  Table  17. 

TABLE  17:  BEARING  FREQUENCIES  OF  MRC310  BEARING 


Bearing  Frequency 

Symbol 

Order 

Train  order  of  rolling 
element 

fr 

0.38 

Relative  rotation  order  of 
rotating  raceway 

ftl 

0.619 

Spin  order  of  rolling 
elements 

fb 

1.98 

Irregularity  order  of 
rolling  element 

fbs 

3.96 

Irregularity  order  of 
rotating  raceway 

fir 

4.95 

Irregularity  order  of 
stationary  raceway 

for 

3.04 

Suppose,  for  example,  that  there  is  a  wear  degradation  in  the  inner  raceway. 
Then  the  presence  of  the  fundamental  f^  tone  with  harmonics  would  be  expected.  Figure 
20(a)  shows  the  spectrum  of  a  good  bearing.  The  6  VdB  fault  is  imposed  at  the  4.95  order 
as  shown  in  Figure  20(b).  Figure  20(c)  shows  the  resulting  increase  in  the  1/1  octave  band 
levels.  The  4.95  order  is  located  at  the  fifth  of  the  1/1  octave  bands.  The  10-5,000  Hz 
broadband  levels  before  and  after  damage  are  almost  unchanged  However,  the  fifth  octave 
band  level  has  a  1 .57  VdB  increase.  This  illustrates  why  a  single  broadband  level  cannot  be 
used  to  determine  the  condition  of  a  machine. 
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Case  MliBearing  Fault 

Type  I  Fire  Pump  (1 1/30/89) 
Pickup  Location  :  MB(FE) 


(a)  Good  Condition 


(b)  Damaged  Condition 
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Figure  20  Simulation  Result  for  Artificial  Bearing  Fault  At  Motor 
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V.  RESULTS  OF  USING  NAVSSES  FREQUENCY  DOMAIN 

DATABASE 

Since  digital  vibration  survey  instruments  are  more  popular  and  easy  to  use  than  other 
types,  only  frequency  spectrum  reading  were  recorded  and  stored  into  the  database.  Thus, 
our  1/1  octave  band  analysis  was  performed  using  frequency  domain  data.  Data  measured 
for  the  Type  I  fire  pumps  in  six  survey  dates  are  available  in  ASCII  format.  The  total 
number  of  fire  pumps  analyzed  was  twenty-seven.  Table  18  shows  the  available  frequency 
data  and  corresponding  survey  date.  Each  set  of  data  contains  a  10  order  and  a  100  order 
spectra.  The  resolution  for  each  set  of  data  is  400  lines.  The  averaged  spectra  using  these 
data  are  shown  in  Figure  21  and  Figure  22  for  each  order  range.  The  thin  lined  curve 
represents  the  mean  levels  and  the  thick  lined  curve  represents  for  the  mean  plus  one 
standard  deviation  levels.  The  10  order  spectrum  was  used  to  compute  the  first  six  1/1 
octave  band  levels  and  the  100  orders  spectrum  was  used  to  compute  the  seventh  and  above 
band  levels.  A  summary  of  1/1  octave  band  analysis  results  for  each  pickup  location  using 
all  the  fire  pumps  is  shown  in  Figure  23.  As  can  be  seen  from  Figure  23,  the  higher 
frequency  levels  (seventh  band  to  ninth  band)  are  about  95  VdB  for  the  motor  bearing  and 
110  VdB  for  the  pump  bearing.  However,  these  higher  frequency  levels  in  Figure  14  using 
time  waveform  tape  data  are  about  80  VdB  for  the  motor  bearing  and  90  VdB  for  the  pump 
bearing.  As  we  reviewed  the  narrow  band  plot  for  the  NAVSSES  frequency  domain  data, 
we  found  that  the  levels  below  54  VdB  have  been  cut  off  and  forced  to  54  VdB.  This  is  why 
the  higher  frequency  levels  using  the  NAVSSES  frequency  domain  data  are  higher  than  the 
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levels  from  using  time  waveform  tape  data.  Figure  24  shows  a  good  example  of  a  cut  off 
spectrum  for  a  motor  bearing  using  the  NAVSSES  frequency  domain  data. 


Table  18:  DATA  COLLECTION  SHEET  FOR  AVAILABLE  FREQUENCY 
DOMAIN  DATA  FROM  NAVSSES  DATABASE 


Survey  Date 

Number  of  Fire 
Pumps 

Number  of  Pickup 
Location 

Number  of  Data 
Set 

12-MAY-93 

6 

5 

30 

23-FEB-92 

5 

5 

25 

21-OCT-90 

5 

5 

25 

07-JUN-90 

4 

5 

20 

26-FEB-90 

1 

5 

5 

30-OCT-89 

6 

5 

30 
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Narrow  Band  Average  for  0-10  Order 

Using  NAVSSES  Database 

Type  I  Fire  Pump 

Each  Position  Use  27  Pumps  for  Averaging 


(a)  Pickup  Location  :  MB(FE) 
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(b)  Pickup  Location  :  MB(CE) 
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Figure  21  The  Averaged  Spectra  for  Type  I  Fire  Pumps  (0-10  Order)  Using 

Frequency  Domain  Data 
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(c)  Pickup  Location  :  PB(CE) 
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(d)  Pickup  Location  :  PB(FE) 
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(d)  Pickup  Location  :  PB(FE/A) 


Figure  21  The  Averaged  Spectra  for  Type  I  Fire  Pumps  (0-10  Order)  Using 

Frequency  Domain  Data  (Cont.) 
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Narrow  Band  Average  for  0-100  Order 

Using  NAVSSES  Database 

Type  I  Fire  Pump 

Each  Position  Use  27  Pumps  for  Averaging 


(a)  Pickup  Location  :  MB(FE) 


40         50         60 
Frequency  (Orders) 


90  100 


i  r 


(b)  Pickup  Location  :  MB(CE) 


i      r 


-i r 


t       r 


10     20     30     40     50     60     70     80     90    100 
Frequency  (Orders) 


Figure  22  The  Averaged  Spectra  for  Type  I  Fire  Pumps  (0-100  Order)  Using 

Frequency  Domain  Data 
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(c)  Pickup  Location  :  PB(CE) 

i            i            i            i            i            i            i 

i             i 

— 

o  120 

LL+O  : 

" 

w 

1 

ft     A 

cp  100 

WL 

9 

Hw>A/u\ 

▼— 

c  ^^»« 

II 

^wV*\^W               iA                                        ^— M..*  . 

m     80 

^-^w^XA^^^^v^-                 ^aa^^^^^W 

TJ 

^~\\aV     S*\^    ^^Sm...       hy>T^        j»aAW  ."" 

•^v^— \^*-*»>^ 

*s 

> 

^-^ — ^~^^ ^v^^          ^■"^W 

■^■^v    _ 

*^  ^^** 

o. 

\y*^r 

>— — N^^—^^/S- 

«w^"^~- 

cd     60 

- 

■D 

> 

1                  1                  1                  1                  1                  1                  1 

1                                      1 

40 

C 

)            10           20           30           40           50           60          70 

Frequency  (Orders) 

30           90 

100 

(d)  Pickup  Location  :  PB(FE) 

i            i            i            i            i            i            i 

i               i 

— 

f  120 

u.+o  : 

— " 

W 

i 

1 

i 

cp  100 

uA 

® 

fL 

T— 

"^    sl\w 

II 

\      A—  ^^^.       A^              ^^.                                                                                                                                                                                                                                                                                A           ■ 

cd     80 

v\^N\AA*-1>a  -k^*<j^   -o^-^-^A 

V^         «*, 

TJ 

S^~/N^^J*V^^V«^/^A^V                 ...        »yv      A 

1^\iMv/* 

W^»-J 

> 

xy^^-^^-/^>^-«s^^>^^^^vvy^>-v_^^^A/vv-'^v^A 

*^»^  B    —  ^      ffm 

N-*^_> 

o^ 

^^ 

m     60 

■ 

TJ 

> 

A  f\ 

i                     i                     i                     i                     i                     i                     i 

1                  1 

40 

( 

)             10           20           30           40           50           60           70 

Frequency  (Orders) 

30           90 

100 

(d)  Pickup  Location  :  PB(FE/A) 

i             i             i             i             i             i             i 

■               i 

— 

g  120 

^l-KJ  : 

" 

lA 

cp  100 

ftw?U 

e 

™     ^^^V           .A                        TV 

Y\ 

nv^A     ^Av 

S     80 

\^>^yV-«v-.           ^w^v^a 

/c*A~-^"~ *^ 

-i  I,   7 

"O 

^-^^-/^>^          V^7^W^-^W»V>*^'            ' **^ 

?v^^v^  ^— v. 

> 

^VaCC-^^^ 

^^V>^^^ 

1 — - 

o 

m    60 

- 

T3 

> 

1                    1                    1                    1                    1                   1                   1 

i                i 

40 

( 

)            10           20           30           40           50           60          70 

80           90 

100 

Frequency  (Orders) 

Figure  22  The  Averaged  Spectra  for  Type  I  Fire  Pumps  (0-100  Order)  Using 

Frequency  Domain  Data  (Cont.) 
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1/1  Octave  Band  Analysis 

Using  NAVSSES  Database 

Type  I  Fire  Pump 

Each  Position  Use  27  Pumps 
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Figure  23  The  Summarized  1/1  Octave  Band  Levels  for  Type  I  Fire  Pumps  Using 

NAVSSES  Frequency  Domain  Data 
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(c)  Pickup  Location  :  PB(CE) 


100  order  B.B.  u+o  Level=1 20.92  dB 
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Figure  23  The  Summarized  1/1  Octave  Band  Levels  for  Type  I  Fire  Pumps  Using 
NAVSSES  Frequency  Domain  Data  (Cont.) 
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Figure  24  The  Narrow  Band  Spectra  for  Type  I  Fire  Pumps  Using  NAVSSES 
Frequency  Domain  Data  Measured  at  MB(FE)  and  MB(CE) 


61 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

Both  time  and  frequency  domain  analyses  were  performed  using  fire  pump  vibration 
data  to  determine  an  appropriate  alert  level.  Cross  Peak  Analysis  in  the  time  domain  and  1/ 
1  Octave  Band  Analysis  in  the  frequency  domain  were  used  to  set  the  alert  threshold 
level(s)  The  following  conclusions  are  drawn: 

Time  Domain  Analysis  (Cross  Peak  Analysis) 

a.  The  measured  peak  envelop  data  at  five  different  pickup  locations  follows  a 
Gaussian  probability  distribution  in  the  VdB  domain  well. 

b.  The  u+2o  value  computed  using  a  peak  envelop  probability  density  function  in  the 
VdB  domain  gives  a  broadband  peak  amplitude  alert  level  rather  than  the  energy  content  in 
the  vibration  signal. 

c.  We  used  12  measured  data  sets  (6  fire  pumps  and  2  different  measurement  dates) 
which  were  available  on  cassette  data  tapes.  To  improve  the  quality  of  averaging,  more 
measured  data  sets  would  be  required 

d.  NAVSSES  does  not  save  time  domain  data  on  any  storage  media  at  this  time. 
However,  for  some  shipboard  machinery,  such  as  low  and  high  pressure  air  compressors, 
time  domain  data  is  required  to  cast  measured  data  in  the  time-frequency  domain  to  detect 
the  possible  faults. 
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Frequency  Domain  Analysis  (1/1  Octave  Band  Analysis) 

a  1/1  Octave  Band  Analysis  (OBA)  uses  nine  frequency  band  alert  levels  which 
provides  more  detailed  information  about  machine  condition  than  the  simple  broadband 
level  measured  by  a  vibrometer.  When  a  significant  change  in  VdB  occurs  in  a  particular 
frequency  band  or  bands,  narrow  band  zoom  mode  analysis  can  be  performed  for  the 
selected  frequency  band(s)  to  identify  the  component(s)  which  may  have  faults. 

b.  OBA  divides  the  frequency  range  into  10  bins  over  10  kHz.  The  VdB  level  in  each 
frequency  bin  is  quite  sensitive  to  changes  in  the  energy  content  of  the  measured  vibration 
signal. 

c.  The  OBA  approach  can  be  easily  incorporated  into  the  Digital  Vibration  Survey 
Instrument  (DVSI). 

d.  Five  sets  (27  fire  pumps)  of  400  line  frequency  domain  DVSI  data  were  received 
from  NAVSSES  and  OBA  was  performed.  The  VdB  levels  in  the  higher  frequency  bins 
were  generally  high.  This  was  caused  by  the  fact  that  the  vibration  signals  which  were 
lower  than  the  VdB  limit(54  VdB)  were  set  to  the  lower  VdB  limit  in  a  specified  dynamic 
range. 

e.  OBA  was  also  performed  using  time  domain  data  obtained  from  cassette  data  tapes. 
The  changes  in  VdB  levels  at  a  particular  frequency  and  the  corresponding  sidebands  were 
introduced  in  narrow  band  FFT.  The  results  are  transferred  into  the  octave  band  domain  and 
the  sensitivity  in  VdB  levels  in  frequency  bins  was  analyzed.  The  results  are  quite 
promising. 

f.  NAVSSES  maintains  400  line  frequency  domain  data  for  all  machineries  in  a 
graphics  postscript  format,  not  in  ASCII  or  digitized  format  It  was  not  possible  to  convert 
this  graphics  format  data  to  ASCII  format. 
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Based  on  current  studies,  the  following  recommendations  are  made: 

a.  In  the  current  machinery  monitoring  practice,  frequency  domain  data  is  commonly 
saved  using  DVSI.  However,  for  some  shipboard  machinery  such  as  low  and  high  pressure 
air  compressors,  time  domain  data  is  required  to  cast  the  measured  data  into  the  time- 
frequency  domain  to  detect  the  possible  faults.  We  suggest  saving  time  domain  data  for  the 
limited  number  of  shipboard  machinery. 

b.  The  1/1  Octave  Band  Analysis  Method  is  suggested  to  implement  in  DVSI. 

c.  NAVSSES  maintains  400  line  frequency  domain  data  for  all  machinery  in  graphics 
postscript  format,  not  in  ASCII  or  digitized  format.  It  was  not  possible  to  convert  this 
graphics  format  data  to  ASCII  format.  We  suggest  saving  future  files  in  ASCII  format  to 
permit  further  analysis  of  the  data. 

d.  It  is  suggested  to  saving/storing  DVSI  data  "as  is",  without  alteration. 
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APPENDIX  A.  MATLAB  PROGRAM  CODE 

A.      ON-LINE  HELP  DOCUMENTATION 

This  section  of  the  guide  contains  all  of  the  MATLAB  functions  that  we  have 

developed  during  this  research.  Each  function  contains  a  purpose,  synopsis  and  description 
in  alphabetical  order.  The  MATLAB  built-in  functions  used  in  these  programs  are 
described  in  the  MATLAB  reference  guide. 

The  code  was  written  for  the  MATLAB  4.0  version.  For  those  who  want  to  use  these 
programs  in  older  versions  need  to  consult  the  MATLAB  reference  guide  for  necessary 
modifications. 
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1.     kurt 
Purpose 

kurtosis  or  fourth  moment  of  vectors  and  matrices. 
Synopsis 

kurtosis=kurt(x) 
Description 

krut  calculates  the  kurtosis  or  fouth  moment  value. 

For  vectors,  kurt(x)  is  the  kurtosis  value  of  the  elements  in  vector  x.  For 

matrices,  kurt(x)  is  a  row  vector  containing  the  kurtosis  value  of  each 

column. 
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2.     logbar 
Purpose 

Plot  1/1  octave  band  bar  plot 

Synopsis 

logbar(Vdb) 
logbar(Vdb,LineStyle) 
logbar(Vdb,LineStyle,LineWidth) 
Description 

logbar(Vdb)  draws  a  1/1  octave  band  bar  plot  using  Vdb  as  the  y  axis 
and  the  ANSI  preferred  center  frequencies  and  pass  bands  as  the  x  axis. 
LineStyle  and  LineWidth  are  optinal  parameters.  The  default  LineStyle 
is  a  solid  line.  See  "plot"  for  more  information  about  LineStyle.  The 
default  LineWidth  is  0.5. 
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3.     mstat 
Purpose 

Statistical  analysis  for  time  domain  data. 
Synopsis 

mstat(y,fs,nbar,id) 

mstat(y,fs,nbar,id,manaly,mid,locat,mregen,pri) 
Description 

mstat(y,fs,nbar,mid,locat,manaly,mthre,z,mregen,pri)  will  analyze  the 

original  data_y  using  the  manaly  method  and  produce  a  statistical 

report  with  a  nbar  pdf  plot. 

mstat(y)  allows  the  user  to  execute  the  program  in  an  interactive  way. 


y 

fs 

nbar 

mid 

locat 

manaly 


mthre 


z 
mregen 


pn 


amplitude  (calibrated  velocity  data) 

sampling  frequency  (Hz) 

number  of  bars  in  Probability.  Density  Distribution  Plot 

type  of  machinery 

pickup  location 

method  of  analysis 

=  1  means  All  data  points  analysis 

=2  means  Overall  peak  points  analysis 

=3  means  Cross  peak  points  analysis 

with  threshold  or  not 

=1  without  threshold 

=2  with  threshold 

percentage  of  threshold 

transform  to  Vdb  or  not 

=  1  without  transform 

=2  transform  to  Vdb  domain 

print  or  not  ('y'es  or  'n'o) 
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4.     npc 

Purpose 

Cross  peak  points. 
Synopsis 

ypc=npc(y) 

Description 

ycp=npc(y)  will  find  the  peak  points  between  each  two  zero-crossing. 
y  original  data 

yep  cross  peak  data 
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5.     npeak 

Purpose 

Overall  peak  points 
Synopsis 

yp=npeak(y) 
Description 

yp=mpeak(y)  will  find  the  overall  peak  points. 

y  original  data 

yp  overall  peak  data 
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6.     nthre 
Purpose 

Threshold 
Synopsis 

yt=nthre(y,z) 
Description 


yt=nthre(y,z)  will  remove  the  points  with  z%  threshold  . 

y  original  data 

yt  data  after  threshold 

z  percentage  of  the  threshold 

If  z=l  then  the  program  will  remove  the  peak  data  small  than 

1%  of  the  maxinum  peak  value. 
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7.     octave 

Purpose 

1/1  octave  band  level. 
Synopsis 

oct=octave(y) 

oct=octave(y,ij,i2) 
Description 

oct=OCTAVE(y)  returns  the  1/1  octave  band  Vdb  values  from  band  1 

to  band  10. 

oct=OCTAVE(y,ii,i2)  returns  the  1/1  octave  band  Vdb  values  from 

band  ij  to  band  i2  and  the  others  bands  will  be  assigned  a  value  of  0. 

y(:,l) :  Vdb  values 

y(:,2) :  frequencies 
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8.     psd 

Purpose 

To  evaluate  the  averaging  VdB  values  using  Hanning  window. 
Synopsis 

VdB=psd(y,nblock) 
Description 

V dB=psd(y, nblock)  returns  the  average  VdB  values  using  time  domain 

data  y.   The  number  of  averagings   is   specified   by   nblock.   The 

recommend  number  of  averagings  is  20.  The  data  matrix  >>  must  be  a  2n 
by  2  matrix. y(:,  I)  contains  time  series  andy(:,2)  contains  corresponding 
velocity  data.  The  hanning  window  is  used  to  damp  out  the  effects  of  the 

Gibbs  phenomenon  The  reference  OVdB  is  10"8  m/sec. 
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9.     skew 
Purpose 

Skewness  or  third  moment. 
Synopsis 

skewness=skew(x) 
Description 

skew  calculates  the  skewness  or  third  moment  value. 

For  vectors,  skew(x)  is  the  skewness  value  of  the  elements  in  vector  x. 

For  matrices,  skew(x)  is  a  row  vector  containing  the  skewness  value  of 

each  column. 
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10.  transfer 
Purpose 

Plot  time  series  and  frequency  spectrum. 
Synopsis 

[Ya,y,fx]=transfer(yo,t,fo,nblock,sec) 
Description 

transfer(yo,t,fo,nblock,sec)  return  three  vectors,  plot  time  series  and 

frequency  spectrum.  The  return  vectors  are: 

Ya  averaged  fft  vector  of  yo 

y  calibrated  velocity  vector 

fx  frequency  vector 

The  input  parameters  include: 

yo  original  amplitude 

t  time  series  corresponds  to  yo 

fo  operating  frequency  (Hz) 

nblock        number  of  data  blocks 

sec  interested  frequency  sections 

The  spectrum  will  indicate  the  peak  position  assign  in  sec. 

For  example,  input  sec=[5  150],  the  program  will  search  the  peak  point 

between  5  Hz  and  1 50  Hz. 
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B.      PROGRAM  CODE 


1.     kurt.m 


function  kurtosis  =  kurt(y) 

%  KURT  Kurtosis. 

%  For  vectors,  KURT(y)  returns  the  kurtosis. 

%  For  matrices,  KURT(y)  is  a  row  vector  containing 

%  the  kurtosis  of  each  column.  See  also  COV 

%  Chao-Shih  Liu  4-26-93 

[n,m]=size(y); 
ifn==l 

y=y'; 

[n,m]=size(y), 
end 
for  i=l:m 

meany(:,i)=mean(y(:,i)).*ones(n,l); 

stdy  (:  ,i)=std(y(:  ,i))*ones(n,  1 ); 
end 

s=((y -meany ). /stdy ) .  A4 ; 
kurtosis  =  sum(s)/n; 
end 
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2.     logbar.m 

function  logbar(Vdb,LineStyle,Line  Width) 

% 

%  logbar  Plot  1/1  octave  band  bar  plot  using  semilog  in 

%  x  axis(freq.)  and  linear  scale  in  Y  axis(Vdb) 

% 

%  logbar(Vdb)  plots  the  semilog  bar  plot 

%  for  1/1  octave  band  analysis. 

% 

%  logbar(Vdb,LineStyle)  plots  the  semilog  bar  plot 

%  for  1/1  octave  band  analysis  with  various 

%  types.  See  "plot"  for  more  about  LineStyle 

% 

%  Line  Width  is  a  optional  parameter  to  control  the  line 

%  width  of  the  bar.  Default  value  is  0.5 


%  Written  By  Liu,Chao-Shih  7/27/93 

if  nargin<2 

lineStyle='-'; 
end 
if  nargin==  1 

linewidth=0.5; 
end 

% 

%  Define  the  upper  and  lower  limit  of  the  pass  bands  (Hz) 

% 

flow=[11.2  22.4  45  90  180  355  710  1400  2800  5600], 

fhigh=[flow(2:l  0)10000], 


%  Determine  the  coordinates  of  the  bars 

% 

xl=flow(:); 

xh=fhigh(:), 

y=Vdb(:); 

[n,m]=size(y), 

n3=n*3; 

yy=zeros(n3+l,m), 

xx=yy; 

xx(l:3:n3,:)=xl(l:n); 

xx(2:3:n3,:)=xh(l:n); 

xx(3:3:n3,:)=xh(l  :n); 

xx=[xl(l);xx], 

yy(l:3:n3,:)=y; 

yy(2:3:n3,:)=y; 

yy=[0;yy], 

set(semilogx(xx,yy,LineStyle),'linewidth',Line  Width) 

end 
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3.      mstat.m 

function  mstat(y,fs,nbar,mid,locat^nanaly,nithre^jnregen,pri) 

% 

%  mstat(y,fs,nbar,mid)locat,manaly,mthre>z,mregen,pri) 

%  will  analysis  the  original  data  y  and  produce  a 

%  statistical  report  with  a  nbar  pdf  plot. 

% 

%  mstat(y)  will  allow  the  user  to  perform  the  program  in  an 

%  interactive  way. 

% 

%  y  :  amplitude  (after  calibration) 

%  fs  :  sampling  frequency  (Hz) 

%  nbar  :  number  of  bars  in  Prob.  Density  Distribution  Plot 

%  mid  :  type  of  machinery 

%  locat :  pickup  location 

%  manaly:  method  of  analysis 

%  =1  means  All  data  points  analysis 

%  =2  means  Overall  peak  points  analysis 

%  =3  means  Cross  peak  points  analysis 

%  mthre  :  with  threshold  or  not 

%  =1  without  threshold 

%  =2  with  threshold 

%  z  :  percentage  of  threshold 

%  mregen:  transform  to  Vdb  or  not 

%  =1  without  transform 

%  =2  transform  to  Vdb  domain 

%  pri  :  print  or  not 

%  ='y'  means  to  pnnt  out  the  result 

% 

%  Liu,Chao-Shih  6/7/93 

% 

if  nargin=l 

fs=input(' Input  the  sampling  frequency  ?  (Hz)'); 

nbar=input(' Input  the  number  of  bar  used  in  the  PDF  plot  ?'); 

mid=input('  Input  the  type  of  the  machinery  ?',V); 

locat=input(' Input  the  pickup  location  ?',V); 

manaly =menu(' Choose  the  method  of  analysis',... 

'All  data  points  analysis',... 

'Overall  peak  analysis',... 

'Cross  peak  analysis'), 

mthre=menu('Type  of  the  threshold  ?',... 

'Without  threshold',... 

'With  n%  threshold'); 

if  mthre=2 

z=input(' Input  the  %  of  the  threshold '),z=z/ 100; 

end 

mregen=menu(' Transform  the  data  ?',... 

'Original  data',... 

'Vdb  Domain'); 
end 

% 
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%  Method  of  analysis 

% 

if  manaly==l 

yp=y. 

mt='All  data  points' 
else  if  manaly==2 
yp=npeak(y); 
mt='Overall  peak', 
else 

yp=npc(y); 
mt=' Cross  peak'; 
end 

% 

%  Type  of  threshold 
% 

if  mthre==l 

tht='  without  threshold'; 
else 

tht=['  with  '  num2str(z*100)  '%  threshold']; 
yp=nthre(yp,z); 
end 
pnumber=length(yp); 

% 

%  Compute  the  first  four  mements  for  linear  domain 

% 

yp=abs(yp); 

mepx=['m='  num2str(mean(yp))]; 

stdpx=['s='  num2str(std(yp))]; 

skewpx=['Skewness='  num2str(skew(yp))]; 

kurtpx=['Kurtosis='  num2str(kurt(yp))]; 

% 

%  Compute  the  first  four  moments  for  Vdb  Domain 

% 

if  mregen==2 

yr=20*logl0(yp/le-8); 

meanyr=mean(yr);stdyr=std(yr); 

skewyr=skew(yr);kurtyr=kurt(yr); 

merpx=['m  =  '  num2str(mean(yr))]; 

stdrpx=['s  =  '  num2str(std(yr))]; 

skewrpx=[ '  Skewness= '  num2str(skew(y  r))] ; 

kurtrpx=['Kurtosis='  num2str(kurt(yr))]; 
end 

% 

%  Plot  the  overall  envelope  PDF 

% 

clg, 

if  mregen==l 

subplot(312) 

[av,au]=hist(yp,nbar);bar(au,av); 

ylabel('Probability  ');xlabel(' Velocity  (m/sec)'); 

title('Original  Data(x)  Envelope'); 
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subplot(3 1 3),axis('off '); 

set(text(0. 3, 1  ,mepx),  'fontn ',  'symbol' ,  'fonta ',  'italic'); 
set(text(0. 3,0. 8,stdpx),'fontn',  'symbol',  'fonta', 'italic'); 
set(text(0.3,0.6,skewpx),  'fontn',  'times',  'fonta',  'italic'); 
set(text(0. 3,0. 4  Jcurtpx), 'fontn',  'times',  'fonta',  'italic'); 
else 

subplot(323) 

[av,au]=hist(yp,nbar);bar(au,av); 
ylabel('Probability'),xlabel(' Velocity  (m/sec)'); 
title('(a)  Linear  Domain'); 
subplot(324) 

|av,au]=hist(yr,nbar),bar(au,av); 
ylabel('Probability');xlabel('Vdb  (0Vdb=le-8  m/sec)'); 
title('(b)  Vdb  Domain'); 
subplot(3 1 3);axis('off '); 

set(text(0. 1,1,  mepx), 'fontn',  'symbol',  'fonta',  'italic'); 
set(text(0. 1,0.8, stdpx), 'fontn', 'symbol', 'fonta', 'italic'); 
set(text(0.1,0.6,skewpx), 'fontn', 'times', 'fonta', 'italic'); 
set(text(0.1,0.4,kurtpx),'fontn', 'times', 'fonta', 'italic'); 
set(text(0. 7,  l,merpx), 'fontn',  'symbol',  'fonta',  'italic'); 
set(text(0. 7,0. 8,stdYpx), 'fontn  Ysymbol'.'fontaYitalic'); 
set(text(0.7,0.6,skewrpx), 'fontn', 'times', 'fonta', 'italic'); 
set(text(0. 7,0. 4,kurtrpx), 'fontn', 'times', 'fonta', 'italic'); 
end 

% 

%  Statistical  analysis  results 

% 

subplot(3 1 1  ),axis('off '); 

tal=[' Statistical  Analysis  Results  for  '  mid], 

set(text(0.5,l,tal),'fontsi',14,'fontw','bold','hori','c'); 

loc=[' Pickup  Location  :  '  locat], 

set(text(0.5,0.8>loc),'fontsi',14,'fontw','bold','hon','c,)> 

ta2=['Data  analysis  by  using  :  '  mt  tht]; 

set(text(0.5,0.4,ta2),'hon','c'); 

block=[' Number  of  data  points  :  '  num2str(pnumber)]; 

set(text(0.5,0.2,block),'hori','c'); 

frange=['The  sampling  frequency  :  '  num2str(fs) '  Hz']; 

set(text(0.5,0,frange),'hori','c'); 

% 

%  Print  the  graph  or  not 

% 

if  nargin<4 

pri=input(Tnnt  the  plot  ?  (y/n)','s'); 
end 
if  pri— 'y' 

set(l,'PaperPosition',[1.25  2  6  8.5]) 

print  -f  1 
end 
end 
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4.     npc.m 


function  ycp=npc(y) 

% 

%  ycp=npc(y)  will  find  the  peak  points  between  each  two 

%  zero-crossing. 

% 

%  y  :  original  data 

%  yep  :  cross  peak  data 

% 

% 

%  Liu,Chao-Shih  4/26/93 

% 

y=yC); 

k=length(y), 

y=y-mean(y)*ones(k,l); 

npc=0;l=l; 

ycp=zeros(k,l); 

for  n=2:k 

if  sign(y(n)*y(n-l))==-l 

h=n-l;npc=npc+l; 

[yp,np]=max(abs(y(kh))); 

np=np+l-l; 

ycp(npc)=y(np); 

l=h+l; 

end 
end 

ycp=ycp(l:npc); 
end 
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5.     npeak.m 


function  yp=mpeak(y) 

% 

%  yp=mpeak(y)  will  find  the  overall  peak  points. 

% 

%  y  :  original  data 

%  yp  :  overall  peak  data 

% 

% 

%  Liu.Chao-Shih  4/14/93 

% 

y=yC); 

k=length(y); 

y=y-mean(y)*ones(length(k),l ); 

yp=zeros(k,l); 

yp(l)=y(l); 

pc=l; 

for  n=2:k-l 

nsl=n-l, 

nl=sign((y(n+l)-y(n))*(y(n)-y(nsl))); 

if  nl==-l  |y(n)=y(nsl) 

pc=pc+l;yp(pc)=y(n); 

end 
end 

yp=yp(l:pc,l); 
end 
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6.     nthre.m 

function  yt=nthre(y,z) 

% 

%  yt=nthre(y,z)  will  remove  the  points  with  z%  threshold 

% 

%  y  :  original  data 

%  yt :  data  after  threshold 

%  z  :  percentage  of  the  threshold 

%  If  z=l  then  the  program  will  not  return  the  peak 

%  peak  data  small  than  1%  of  the  maxinum  peak  value. 


%  Liu,Chao-Shih  4/14/93 

% 

y=sort(y(:)); 

ly=length(y); 

maxy=max(y); 

miny =z*0.0 1  *maxy; 

forl=l:ly 

if  yO^miny 

po=l; 

yt=y(po:ly,l); 

break 

end 
end 
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7.     octave.m 

function  oct=octave(y,i  1  ,i2) 

% 

%  OCTAVE  Octave  band 

% 

%  oct=OCTAVE(y)  returns  the  1/1  octave  band  Vdb  values  from 

%band  1  to  band  10. 

% 

%  oct=OCTAVE(y,il,i2)  returns  the  1/1  octave  band  Vdb  values 

%  from  band  i  1  to  band  i2  and  the  others 

%  bands  will  be  assigned  to  0. 

% 

%y(:,l):  Vdb  values 

%  y(:,2) :  frequencies 

% 

%Chao-Shm  Liu  7/13/93 

% 

fc=[16  31.5  63  125  250  500  1000  2000  4000  8000]; 
flow=[11.2  22.4  45  90  180  355  710  1400  2800  5600], 
fhigh=[22.4  45  90  180  355  710  1400  2800  5600  11200], 
np=max(size(y)); 
oct=zeros(10, 1); 

% 

%  Default  the  ll  and  i2 
% 

if  nargin==  1 
il  =  l;i2=10, 
end 

% 

%  if  Vdb  lower  than  0  db,  force  the  value=0db 

% 

for  k=  1  :np 

ify(k,l)<0 

y(k,l)=0, 

end 
end 

% 

%  convert  to  linear  velocity  value 

% 

vref=le-8; 

y(:,l)=10A(y(:,l)/20)*vref; 

fmax=y(np,2);df=fmax/(np- 1 ), 

% 

%  Calculate  the  1/1  octave  band  levels  in  linear  domam 

% 

df2=df/2, 

fupper=y(:,2)+df2,  %  Upper  limit  of  frequency 

flower=y(:,2)-df2,  %  Lower  limit  of  frequency 
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sum=zeros(i2-il,l); 
fork=il:i2; 

mark  1  =fupper-flow(k); 
mark2=flower-flow(k), 
mark3=fupper-fhigh(k); 
mark4=flower-fhigh(k); 
for  j=l:np 

if  mark4(j)>=0  %  larger  than  the  band(k) 
break 
else 
if  markl(j)<=0  %  smaller  than  the  band(k) 

%  do  nothing 
else  %  cover  the  band(k) 
ifmark3(j)>0 

if  mark2(j)<0  %  band(k)  smaller  than  df 
sum(k)=y(j,l)A2*(fhigh(k)-flow(k)); 
break; 

else  %  partial  band  in  band(k+l) 
sum(k)=sum(k)+y(j,l)A2*abs(mark4(j)X 
end 
else 

if  mark2(j)<0  %  partial  band  in  band(k-l) 
sum(k)=y(j.  1  )A2*markl  (j); 
else 

sum(k)=sum(k)+y  (j,  1  )A2*df, 
end 
end 
end 
end 
end 

%  Convert  the  octave  band  levels  from  linear  scale  to  VdB  levels 
oct(k)=20*log  1 0(sqrt(sum(k))/vref); 
end 

%  Plot  the  results 
bar(oct) 

xlabel('Number  of  Center  Frequency ') 
ylabel('Vdb  value  (0Vdb=le-8  m/sec)') 
axis([il-l  i2+l  40  120]) 
fork=U:i2 

octext=spnntf('%5.  If  \oct(k)), 

set(text(k,oct(k)+2,octext),,fontw^'bold^'ho^2^, center') 
end 
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8.     psd.m 


function  Vdb=psd(y,nblock) 

%  Vdb=PSD(y,nblock) 

% 

%  This  function  will  compute  the  average  Vdb  value  by 

%  using  hanning  to  nblock  data  set. 

% 

%  y  :  original  data  matrix  (time.velocity) 

%  nblock  :  number  of  data  blocks 

%  Ya  :  average  fft  Vdb  vector  of  y(0Vdb=le-8  m/sec) 

%  fx  :  frequency  vector 

% 

%  Liu,Chao-Shih  12:11PM  7/14/93 

% 

%  standardize  the  data  matrix 

n=size(y,2); 

ifn~=2 

y=y'; 

end 

t=y(:,i);y=yC,2); 

dt=t(2)-t(l); 

fs=l/dt; 

vref=le-8; 

% 

%  Compute  the  fft  data 

% 

n=length(y)/nblock,n2=n/2; 

nl=l; 

w=hanning(n); 

for  nb=  1 :  nblock 

nh=nb*n; 

Y(l  :n,nb)=fft(w  .*y(nl:nh,l))*dt; 

nl=nh+l; 
end 
if  nblock~=l 

Y=(sum(abs(Y)')/nblock),> 
else 

Y=abs(Y); 
end 

Ya=20*loglO(2*Y/vreO; 
f=fs/2*(0:n2)/n2; 
Ya(n2+2:n)=[]; 
fx=f(2:n2+l)'; 
df=fx(2)-fx(l); 
fx=[0;fx], 
Vdb=[Ya,fx], 
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skew.m 

function  skewness  =  skew(y) 

%SKEW  Skewness.  For  vectors,  SKEW(y)  returns  the  skewness. 

%  For  matrices,  STD(y)  is  a  row  vector  containing 

%  the  skewness  of  each  column.  See  also  COV. 

%  Chao-Shih  Liu  4-26-93 

[n,m]=size(y); 
if  n==l 

y=y'; 

[n,m]=size(y); 
end 
for  i=l:m 

meany(:,i)=mean(y(:,i)).*ones(n,l); 

stdy(:,i)=std(y(:,i))*ones(n,  1 ); 
end 

s=((y-meany)./stdy).A3; 
skewness  =  sum(s)/n; 
end 
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10.   transfer.m 

function  [Ya>y,fx]=transfer(yo,t,fo,nblock,sec) 

%  [Ya,y)fx]=transfer(yo,t>fo,nblock,sec) 

% 

%  This  function  will  plot  time  &  frequency  domain 

% 

%  yo  :  original  velocity  vector  (before  calibration) 

%  t :  time  vector 

%  fs  :  sample  frequency 

%  fo  :  operation  frequency  =  RPM/60 

%  nblock  :  number  of  data  blocks 

%  Ya  :  average  fft  vector  of  y 

%  fx  :  frequency  vector 

%  y  :  original  velocity  vector  (after  calibration) 

% 

%  Lru.Chao-Shih  1:58  7/10/93 

% 

dt=t(2)-t(l), 

fs=l/dt; 

itext=input(' Input  the  title  of  this  record  =>','s'); 

locat=input(Tnput  the  location  of  data  collect'.'s'); 

y=yo; 

scale=menu(' Choose  the  method  of  X  scale' /Frequency  (Hz) '.'Order'); 

vref=le-8; 


%  Standanzed  the  size  of  y 

% 

[xd,yd]=size(y); 

if  yd~=l 

y=y';t=t'; 

end 
ny=length(y); 

% 

%  Compute  the  fft  data 

% 

n=ny/nblock;n2=n/2; 

nl=l; 

win=menu('Type  of  Window  System', 'Without  window  system' 

'Harming  window', 'Hamming  window'), 
if  win=l 

w=ones(n,l); 
elseif  win=2 

w=hanning(n); 
else 

w=hamming(n), 
end 
for  nb=l:  nblock 

nh=nb*n; 

Y(  1  :n,nb)=fft(w.  *y  (nl  :nh,  1  ))*dt ; 

nl=nh+l; 
end 
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if  nblock~=l 
Y=sum(abs(Y)')/nblock; 
else 

Y=abs(Y); 
end 

Ya=20*loglO(2*Y/vref), 
f=fs/2*(0:n2)/n2; 
Ya(n2+2:n)=[]; 
if  scale=l 
fx=f(2:n2+l); 
df=fx(2)-fx(l); 

xt= 'Frequency  (Hz)';unit='  Hz  '; 
else 

fx=f(2:n2+l)/fo; 
df=fx(2)-fx(l); 

low=ceil(0.75/df),high=ceil(  1 .25/df); 
[  Ymax,Ypos]=max(Ya(low:high));fscale=fx(low+Ypos- 1 ), 
fx=fx/fscale;fo=fo*f scale; 

xt=['Frequency  (Normalized  Orders) '];unit='  Order  '; 
end 

df=fx(2)-fx(l); 
sec=ceil(sec/df); 
ifsec(l)==0 
sec(l)=l; 
end 

for  mk=  1 :  length(sec)- 1 

[mav(mk),pos(mk)]=max(Ya(sec(mk):sec(mk+l))); 
pos(mk)=pos(mk)+sec(mk); 
end 


%  Plot  time  domain  and  frequency  domain 

% 

figure(  1  );clf  ;subplot(2 1 2); 

tf=1024; 

ifny<1024 

tf=ny; 
end 

plot(t(l:tf,l))y(l:tf,l),'r,);axis([t(l,l)t(tf,l)min(y)max(y)]); 
titleCTimeDomain');xlabel('Time  (sec)');ylabel('Velocity(m/sec)');grid 
% 

%  General  information 
% 

subplot(2 1 1  );axis('off ') 
mfo=['The  operation  condition  =  '  num2str(fo*60)... 

'  RPM  (=  '  num2str(fo) '  Hz)']; 
mloc=['Data  pickup  location  =  '  locat]; 
mfs=['Data  sampling  frequency  =  '  num2str(fs) '  Hz']; 
mdp=[Total  number  of  sampling  data  points  =  '  num2str(length(y))]; 
set(text(0.5,l,itext),'fontsi',15,'fontw', 'bold', 'horiz', 'center') 
text(0.05,0.9,mloc);text(0.05,0.82,mfo); 
text(0.05,0.74,mfs);text(0.05,0.66,mdp); 

text(0,0.5,'*  The  time  domain'); 

mt=['Mean  value  =  '  num2str(mean(y)) '  m/sec']; 
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text(0.05,0.4,mt); 

st=[' Standard  deviation  =  '  num2str(std(y)) '  m/sec']; 

text(0.05,0.32,st); 

text(0. 05,0.24, ['The  plot  only  shown  '  num2str(tf) '  data  points']), 

pause 

figure(2),clf  ;subplot(2 1 2); 
plot(fx,Ya(2:n2+l),'r'); 
upper=(ceil((max(mav))/10)+l  )*  10; 
axis([0  fx(n2)  20  upper]), 

title( 'Frequency  Domain');xlabel(xt),ylabel(' Velocity  in  Vdb');grid 
for  lop=l  :length(sec)-l 
text(fx(pos(lop)-l),mav(lop),num2str(lop)); 
end 


%  General  information 

% 

subplot(211);axis('off) 

mfo=['The  operation  condition  =  '  num2str(fo*60).. 

'  RPM  (=  '  num2str(fo) '  Hz)']; 
mloc=['Data  pickup  location  =  '  locat], 
mfs=['Data  sampling  frequency  =  '  num2str(fs) '  Hz']; 
mbk=['Number  of  data  block  averages  =  '  num2str(nblock)]; 
mdp=['Block  size  =  '  num2str(length(y)/nblock) '  points/block' 
set(text(0.5,l,itext),'fontsi',15,'fontw','bold','horiz','center'); 
text(0.05,0.9,mloc),text(0.05,0.82,mfo); 
text(0.05,0.74jnfs),text(0.05,0.66,mbk),text(0.05,0.58,mdp); 
text(0,0.5,'*  The  frequency  domain'); 
for  lop=l  :length(sec)-l 

as=[num2str(lop),'  =  (',num2str(fx(pos(lop)-l )),... 

unit,',  ',num2str(mav(lop)),'  db)']; 

if  lop>4 

xpos=0.55;ypos=0.8-0.08*lop; 

else 

xpos=O.05;ypos=0.48-0.08*lop; 

end 

text(xpos,ypos,as); 
end 

morder=['l  order  =  '  num2sU-(fo) '  Hz']; 
text(0.05,0.084norder), 

mref=['The  reference  0  Vdb  =  '  num2str(vref) '  m/sec']; 
text(0.05,0,mre0; 
% 

%  zooming  up  the  plot  ? 
% 

zoom=  1 , 
while  zoom~^=3 
subplot(212) 
zoom=menu('zoom  up  ?','0-10  order',... 

'Use  mouse  to  select  the  range', 'No'); 

if  zoom==  1 

axis([0  10  20  upper]); 

else  if  zoom==2 

axis([0  fx(n2)  20  upper]); 
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[xs,ys]=ginput(2); 

xrange=sort(xs); 

if  xrange(l)<0 

xrange(l)=0; 

end 

axis([xrange(l)  xrange(2)  20  upper]); 

else 

break 

end 
end 
% 

%  Print  out  the  plot 
% 
uiprint=[  'set(l,"PaperPosition",[0.25  1  8  9.25]);'  'print  -fl;' ... 

'set(2,"PaperPosition",[0.25  1  8  9.25]);'  'print -f2;'  'end']; 
uicontrol('style', 'push', 'units', 'normal', 'pos',[0.9  0.9  0.08  0.06], . 

'string',  'print',  'call  \uiprint); 
uicontrolC style', 'push', 'units', 'normal', 'pos', [0.9  0.82  0.08  0.06], 

'string', 'exit', 'call', 'end'); 
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APPENDIX  B.  EXAMPLES  OF  USING  MATLAB  PROGRAM  IN 

PC486 

MATLAB  allows  the  user  to  write  a  script  file.  A  script  file  is  an  external  file  that 
contains  a  sequence  of  MATLAB  statements.  The  following  example  will  show  the  script 
files  used  to  perform  the  analysis. 

A.      CROSS  PEAK  ANALYSIS  EXAMPLE 

In  this  example,  a  CG-59  fire  pump  measured  at  MB(FE)  data  set  is  used  in  the 
analysis.  The  file  name  of  the  data  file  \s  flpl.asc.  The  MATLAB  command  is  anchored 
by». 

» load  flpl.asc 

»y=flpl(:,2); 

»t=flpl(:,l); 

»  mstat(y) 
The  user  can  then  interact  with  the  program  to  select  the  method  of  analysis.  Instead 

of  this  interactive  process,  the  user  can  use  the  following  to  execute  the  program  in  an 

automatic  way. 

» load  flpl.asc 
»y=flpl(:,2); 

»t=flpl(:,l); 

»  ypc=npc(y); 

»  yt=nthre(ypc,l), 

»  mid='CG-59  Fire  Pump', 

»locat='MB(FE)', 
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»  mstat(yt,  1 0000, 1 0,2,mid,locat,3,2, 1 ,2,'y ') 

B.      1/1  OCTAVE  BAND  ANALYSIS  EXAMPLE 

Here,  proceeding  as  in  the  above  example,  a  CG-59  fire  pump  measured  at  MB(FE) 
data  set  is  used  to  perform  the  1/1  octave  band  analysis.  The  file  name  of  the  data  file  is 
flpl.asc.  The  MATLAB  command  is  anchored  by  ». 

» load  flpl.asc 

»  yh=flpl(:,2);  %  yh  is  a  time  series  sampled  at  10  kHz. 

»  yl=yh(l  :5:length(yh));  %  yl  is  a  time  series  sampled  at  2  kHz. 

»t=flpl(:,l); 

»  VdBh=psd(yh,20); 

»  VdBl=psd(yl,20); 

»  octh=octave(VdBh,7,9); 

»  octl=octave(Vdbl,l,6); 

»  octband=[octh  octl]; 

» logbar(octband); 
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APPENDIX  C.  FIGURES  OF  STATISTICAL  ANALYSIS  RESULTS 

This  section  contains  figures  of  the  statistical  analysis  results  with  and  without 
considering  damaged  data  sets. 
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Statistical  Analysis  Results 

Type  1  Fire  Pump 
Pickup  Location  :  MB(FE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  1.965e+005 
The  sampling  frequency  :  1e+004  Hz 

(a)  Linear  Domain                                 (b)  VdB  Domain 
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Figure  C.l  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
MB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  I  Fire  Pump 
Pickup  Location  :  MB(CE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  2.173e+005 
The  sampling  frequency  :  1e+004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 
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Figure  C.2  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
MB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 


96 


Statistical  Analysis  Results 

Type  I  Fire  Pump 
Pickup  Location  :  PB(CE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  3.129e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.3  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
PB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  I  Fire  Pump 

Pickup  Location  :  PB(FE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  2.614e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.4  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
PB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  I  Fire  Pump 
Pickup  Location  :  PB(FE/A) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  3.105e+005 
The  sampling  frequency  :  1e+004  Hz 
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(b)  VdB  Domain 
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Figure  C.5  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
PB(FE/A)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  II  Fire  Pump 
Pickup  Location  :  MB(FE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  4.733e+004 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.6  The  Statistical  Analysis  Results  for  Type  II  Fire  Pumps  Measured  at 
MB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  II  Fire  Pump 
Pickup  Location  :  MB(CE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  4.835e+004 
The  sampling  frequency  :  1e+004  Hz 
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(b)  VdB  Domain 
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Figure  C.7  The  Statistical  Analysis  Results  for  Type  II  Fire  Pumps  Measured  at 
MB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  II  Fire  Pump 
Pickup  Location  :  MB(CE/A) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  3.089e+004 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.8  The  Statistical  Analysis  Results  for  Type  II  Fire  Pumps  Measured  at 
MB(CE/A)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  II  Fire  Pump 
Pickup  Location  :  PB(CE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  3.241  e+004 
The  sampling  frequency  :  1  e+004  Hz 
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Figure  C.9  The  Statistical  Analysis  Results  for  Type  II  Fire  Pumps  Measured  at 
PB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  II  Fire  Pump 

Pickup  Location  :  PB(FE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  4.056e+004 
The  sampling  frequency  :  1e+004  Hz 
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Figure  CIO  The  Statistical  Analysis  Results  for  Type  II  Fire  Pumps  Measured  at 
PB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  UMB 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 

Number  of  data  points  :  2.1e+005 

The  sampling  frequency  :  1e+004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 


0.3 


.J  0.2 

CO 
JQ 

2 

CL  0.1 


_n 

n~ 

0.005         0.01 
Velocity  (m/sec) 


^1=0.0006716 
0=0.0009258 
Skewness=4.985 
Kurtosis=36.52 


0.015 


50  100  150 

VdB  (0VdB=1e-8  m/sec) 


\i  =  92.\2 
a  =  8.687 
Skewness=0.02892 
Kurtosis=2.953 


Figure  C.ll  The  Statistical  Analysis  Results  for  Type  in  Fire  Pumps  Measured  at 
UMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  LMB 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  2.176e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.12  The  Statistical  Analysis  Results  for  Type  III  Fire  Pumps  Measured  at 
LMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 
Pickup  Location  :  LMB(A) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  2.845e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.13  The  Statistical  Analysis  Results  for  Type  III  Fire  Pumps  Measured  at 
LMB(A)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  UPB 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  1.367e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.14  The  Statistical  Analysis  Results  for  Type  III  Fire  Pumps  Measured  at 
UPB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  LPB 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 

Number  of  data  points  :  1 .99e+005 

The  sampling  frequency  :  1e+004  Hz 
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Figure  C.15  The  Statistical  Analysis  Results  for  Type  III  Fire  Pumps  Measured  at 
LPB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  I  Fire  Pump 

Pickup  Location  :  PB(CE)  w/o  89#6,90#3,90#6 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  2.704e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.16  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  I  Fire 
Pumps  Measured  at  PB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  I  Fire  Pump 

Pickup  Location  :  PB(FE)  w/o  90#3 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  2.489e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.17  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  I  Fire 
Pumps  Measured  at  PB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  II  Fire  Pump 

Pickup  Location  :  MB(CE)  w/o  #1 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  3.238e+004 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.18  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  II  Fire 
Pumps  Measured  at  MB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  UMB  w/o  11/24/#1 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 

Number  of  data  points  :  2.05e+005 

The  sampling  frequency  :  1e+004  Hz 
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Figure  C.19  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  III  Fire 
Pumps  Measured  at  UMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  LMB  w/o  11/24/#1,3/10/#2 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  1 .873e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.20  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  III  Fire 
Pumps  Measured  at  LMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 

Type  III  Fire  Pump 

Pickup  Location  :  UPB  w/o  11/24/#1 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points  :  1.245e+005 
The  sampling  frequency  :  1e+004  Hz 
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Figure  C.21  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  III  Fire 
Pumps  Measured  at  UPB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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